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Amendments to the Specification: 

Page 1, line 10, please amend the following paragraphs: 

COMPUTER PROGRAM LISTING APPENDIX 

The computer program listing appendix attached hereto consists of two (2) identical 
compact disks, copy I and copy 2, each containing a listing of the software code for 
embodiments of components of this invention. Each compact disk contains the following files 
(date and time of creation, s ize in b>tes, filename): 

Directory of D:\ 

05/30/2003 09:09 AM 0711115B,txt 

05/30/2003 09:11 AM W4 071in5C.txt 

05/30/2003 09:12 AM : b024 0711115D,txt 

05/30/2003 09:18 AM 1,361 0711115E.txt 

05/30/2003 09:19 AM 335 0711115Rtxt 

05/30/2003 09:20 AM 649 0711115G.txt 

05/30/2003 09:08 AM : — 3^9S9 071115A.txt 

05/30/2003 09:05 AM 38,253 071119.txt 

8 Fil e (s) 51,167 bytos 

0 Dir(s) 0 byt e s fr ee 

Th e contents of the compact disk are a part of th e pre se nt disclosur e , and ar e incorporat e d 
by r e f e r e nc e h e r e in in th e ir e ntireti es . 
COPYRIGHT NOTICE 

A portion of th e disclosur e of this pat e nt document contains mat e rial that is subj e ct to 
copyright prot e ction. Th e copyright owner has no objection to th e facsimile reproduction by 
anyone of th e pat e nt docum e nt or th e pat e nt disclosure as it app e ars in th e Patent and Trad e mark 
Offic e patent fil e or r e cords, but oth e rwis e r e serv e s all cop>Tight rights whatso e v e r. 

Page 4, line 10, please amend the following paragraph: 

There are a number of causes of audio distortion v^Heh -that involve the structure and 
operation of the voice coil transducer 300. At high signal levels, voice coil transducers become 
very distorting. This distortion is largely caused by the nonlinear it ies in the coil motor factor, in 
the restoring force of the coil/diaphragm assembly suspension, and the impedance of the coil. 
Other nonlinear effects also contribute to the distortion. Nonlinear effects are an intrinsic part of 
the design of voice coil transducers. 
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Page 5, line 7, please amend the following paragraph: 

FIG. 4 represents data for actual large signal (LS) parameters of a transducer from a small 
desktop stereo system, model name: Spin70, manufactured by Labtec. The large signal 
parameters shown in FIG. 4 were obtained using a commercially available laser metrology 
system (Klippel GMBH). The magnitude of Bl is shown by curve 401 as a function of the 
displacement x of the coil/diaphragm assembly from the no-signal equilibrium position, which 
is indicated in FIG. 4 by a zero on the horizontal axis; at that position, no elastic restoring force 
is applied to the -coil/diaphragm assembly. The unit for Bl is Newton I Ampere (or N I A), The 
highly non-constant nature of the Bl factors of commercial voice coil transducers is recognized 
in the current art. As the audio signal increases in magnitude, the coil tends to move away from 
the region of maximal static magnetic field, and the motor factor decreases, thus effecting a less 
uniform coil movement and distorting the sound wave. 

Page 5, line 23, please amend the following paragraph: 

Referring to FIG. 3, as pointed out above, the cone suspension is axially symmetric and 
typically includes two parts: a corrugated suspension near the coil, typically referred to as the 
spider 307, and the surround 308 connecting the large end of cone 302 to the frame 301 of the 
speaker. These two suspensions together act as an effective spring, which provides a restoring 
force to the coil/diaphragm assembly and determines the equilibrium position of the assembly to 
which it relaxes when not being driven. This effective spring restoring force is again a highly 
non-constant function of coil/cone axial position x ; that is to say, the effective spring stiffhess 
varies significantly as a function of x . In FIG. 4 curve 402 shows a plot of K, the spring 
stiffhess, as a function of x for the speaker transducer mentioned above. Spring stiffhess K is 
expressed in units of N Imm (-i.e. Newtons per millimeter). 

Page 6, line 13, please amend the following paragraph: 
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where m is the mass of the assembly plus a correction for the mass of air being moved; 

R^^ represents the effective drag coefficient experienced by the -assembly, mainly due to air 

back pressure and suspension friction; K{x) is the position dependent effective spring stiffness 
due to the elastic suspension; Bl{x) is the position dependent motor factor; and i{t) is the time 
dependent voice-coil current, which responds to the input audio signal and constitutes the control 
variable. These terms are related to the industry standard linear model (small signal) parameters 
- namely, the Thiele-Small parameters, which are as follows: 

Page 7, line 4, please amend the following paragraphs: 

The second order differential equation (1) would be straightforward to solve, but for the 
nonlinearities in the elastic restoring force and in the motor force terms; these nonlinearities stem 
from the x dependence of K{x) and Bl{x)^ and- they preclude a closed-form analytical solution 
in the general case. Although approximations can be made, it is_difFicult to predict the response 
of a system under all conditions, and thus to create a robust control system. 

Further nonlinearities arise due to other electrodynamical effects caused by the 
application of the audio signal to the transducer voice-coil. Typically, current is supplied to the 
coil by converting the audio information into a voltage, V(t), which is imposed across the 
terminals of the voice coil. However, the resulting coil current varies both out of phase and 
nonlinearly with this voltage. The phase lag arises both because the voice coil's effective 
impedance has a reactive component, and because the electromechanical transduction of the coil 
current into coil motion through the -static magnetic field induces a Back-ElectroMotive Force 
(BEMF) voltage term in the coil circuit. 

Page 7, line 30, please amend the following paragraph: 

— Whef=e -where the BEMF is represented by the second term on the left hand side (a product of 
Bl(x) and coil velocity). The Ohmic resistance of the coil is . The coil's effective 

inductance, Z,(x) , is a function of x because it depends upon the instantaneous position of the 
coil relative to the magnetic pole structure and its airgap. In FIG. 4 curve 403 shows a typical 
plot of the position dependence of coil inductance L^ix) at low audio frequencies. The units of 

4 are mH (milli-Henries), and the values of shown in curve 403 have been multiplied by a 
factor 10 to render the graph more readable. 
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Page 9, line 7: equation (4) 




(4) 



m m 



m m 



Page 10, line 10, please amend the following paragraph: 

Accordingly, there are several factors described above that significantly affect the ability 
to provide -accurate sound from a conventional audio reproduction system. Some of the issues 
can be addressed by improving the circuitry through digital means; but even with the digital 
circuitry to handle the signal shaping, the transducer itself has significant nonlinearities that can 
never be addressed adequately by shaping the input signal to the transducer. Therefore, what is 
needed is a system that controls the transducer in such a manner that optimum linear sound is 
provided. Such a system should also be easy to implement, cost effective, and easily adaptable 
to existing systems. The present invention provides a control system for a transducer to provide 
linear sound, and the present invention also provides an integrated audio reproduction system. 

Page 13, line 18, please amend the following paragraph: 

FIG. 10 is a block diagram of the feedback Hnearization process using the control law 
given by equation f^ (40), which provides transduction corrections along with a linear spring 
constant (suspension stiffness) whieh -that is electronically added; 

Page 14, line 2, please amend the following paragraph: 

FIG. 14 illustrates Power Spectrum Distribution simulation curves showing the effect of 
the transduction corrections (spring stiffness and motor factor correction) upon harmonic 
distortion for a single- 100/fe tone input, both with and without BEMF and nonlinear 



inductance in the physical model of the Labtec Spin 70 transducer; 

Page 14, line 1 1, please amend the following paragraph: 

FIG. 16 illustrates simulated waveforms of the coil/diaphragm axial position versus time 
in the presence of a single-tone excitation, both with and without electronically restored effective 
spring stifBiess, showing that without such restoration the cone may drift from its equilibrium 
position and reach its limit of excursion; 

Page 14, line 20, please amend the following paragraph: 
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FIG. 18 is a graph of the simulated phase lag between coil voltage and coil current as a 
function of audio frequency at low frequencies, which is almost entirely due to BEMF; 

Page 15, line 10, please amend the following paragraph: 

FIG. 24 illustrates a plot of suspension stiffness K in Newtons/mm together with a plot 
of Bl in Newtons/amp, both of which are plotted against for the same Labtec Spin 70 
transducer data; 

Page 15, line 16, please amend the following paragraph: 

FIG. 26 shows a curve whieh -that illustrates the variation of 4 with position at 43kHz 
for a Labtec Spin70 transducer; 

Page 15, line 24, please amend the following paragraph: 

FIG. 31 is a block diagram for a circuit which that, together with parameter estimation, 
measures transducer coil inductance via a supersonic probe tone and reference RL circuit; 

Page 15, line 30, please amend the following paragraph: 

FIG. 33 shows a curve w4H€h -that illustrates the variation of C^^^,^,^^ with V^^^, for driven 

measurement; C^;.^^,/,^ measured in arbitrary units obtained using the method described in 

Detail 12; 

Page 16, line 9, please amend the following paragraph: 

FIG. 37 is a cross-sectional view of a speaker transducer ^v^4H€h -that includes an IR-LED 
diode and an associated PIN diode, mounted on the back side of an audio transducer of the type 
shown in FIG. 3, as part of an optical position detection system; 

Page 16, line 27, please amend the following paragraph: 

FIG. 43 shows the details of HW and ISR operations for the S calibration in step 1 1504 
of FIG. U&42; 

Page 17, line 1, please amend the following paragraphs: 

FIG. 46 illustrates the operations of process 1 1203 of FIG. 45 that are spawned as a result 
of enabling sampling clock and ISR in step 1 1202 of FIG. 44345; 

FIG. 47 shows a flow diagram of the ISR 1 1303 of FIG. +4^46; 
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Page 18, line 16, please amend the following paragraph: 

FIG. 66 shows a low-frequency portion of the same power spectrum distribution shown 
in FIG. 65, displaying multiple harmonics of the 60 Hz tone, with spectra depicted both with and 
without correction. 

DETAILED DESCRIPTION OF THE EMBODIMENT(S) 

Page 19, line 1, please amend the following paragraph: 

An enabling invention in the area of control engineering 501 was the linearization 
method for dynamical equations 504 used in modeling physical systems to be controlled, such as 
actuators and transducers. This method relies on finding the control equation for the non-linear 
part of the dynamical equation and substituting this into the fiill equation. The application of this 
method to a second order differential equation 505 shows that a non-linear second order ordinary 
differential equation can be linearized by solving the control equation for the non-linear first 
order differential equation, provided the second order and first order differential terms are linear. 
This is a general method for linearizing such differential equations, and covers the application to 
the control of all actuators and transducer systems that can be modeled in full, or in part, by such 
an equation. The application of the linearizing method 505 to an equation with non linearities 
dependent on one state variable 506 shows that only one state variable is required for 
linearization. The application of 506 relies on positional sensing. That is to say, neither the 
velocity, nor the acceleration, nor the instantaneous driving force state variables are required in 
order to linearize the process. Position dependent sensing and feedback linearization can be used 
with many classes of a-non-linear motors and actuators. 

Page 19, line 25, please amend the following paragraphs: 

In the present work it was discovered that there are multiple processes in a sound 
reproduction system, that each process can influence the performance of other processes, that 
each process has non-linearities that must be considered in the design of a control paradigm loop. 
and that each control paradigm loop must have a sufficient number of state measurements which 
must be measured with sufficient discrimination against noise and with sufficient speed to 
control the process. It was further discov e red that control of on e proc e ss must b e stabl e in th e 
presenc e of oth e r proc e ss e s. 

Control of multiple processes with multiple control paradigms loops can be aff e ct e d 
effected if the criteria for sufficiency is met for each control paradigm loop. It has been 



Attorney Docket: 3141P 

discovered that for the correction of non-linear transduction a necessary condition for control is a 
positional state measurement, in distinction to the motional measurements of prior art. The 
positional state measurement must be of sufficiently low- noise and latency and of sufficiently 
high speed, or bandwidth^ to effect the control while not adding unacceptable noise to, nor 
engendering instability in, the sound output. Multiple positional measurements can be used to 
estimate the positional state for the purpose of transducer linearization. 

Page 20, line 30, please amend the following paragraph: 

The method and system comprise providing a model of at least a portion of the audio 
transducer system and utilizing a control engineering technique in the time domain to control an 
output of the audio transducer system based upon the model. In the present invention a method to 
determine, in real-_time, the nonlinear parameters of the transducer from measurement of internal 
state parameters of the transducer is provided. In particular the electrical properties of the voice 
coil can be used as a measure of positional state and a predictor of the major non-linearities of 
the transducer. "Real-Jime" in this context means with sufficiently bandwidth low latency to 
effect control. 

Page 21, line 12, please amend the following paragraphs: 

It has been discovered that in an audio reproduction system, the overall process of 
converting audio information into sound can be considered as consisting of three processes. First, 
conditioning of the audio signal to produce the transducer drive signal; second, the transduction 
of the drive signal into a diaphragm motion moving an air mass; and third, the conditioning of 
the moving air mass to provide an output sound. Thus, an audio transducer can be defined as: 
signal conditioning/transduction/sound conditioning. FIG. 6 illustrates a block diagram of audio 
reproduction system 11 00 in accordance with these processes. As is seen, a signal conditioning 
process 1 102 takes an audio signal 1101 (digital or analog) and performs signal conversion, 
amplification, filtering and frequency partitioning to provide a drive signal 1 103. The drive 
signal 1 103 is provided to a transduction process 1 104. The transduction process 1 104 typically 
utilizes a plurality of transducers, and results in diaphragm motion 1 105^ which drives an air 
load. A sound conditioning process 1 106, which may include effects from a speaker enclosure 
and an extended audio environment, acts on the air load driven by the diaphragm motion 1 105 to 
provide the perceived sound 1 107. 

Distorting factors due to nonlinear effects influence all of these processes. These factors 
arise in the relationship between the audio signal as a voltage and the drive current in the coil 
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(transconductance), and in the electro-magneto-mechanical (henceforth abbreviated 
"electromechanical") effects involving the moving-coil motor. Nonlinear effects resulting from 
sound conditioning are much smaller in normal operating conditions, and are thus neglected in 
the physical model described in this section, and in the control model based upon it and 
described in Detail 2. But these nonlinear acoustical effects, along with other higher-order effects 
described and then neglected in this section, can in principle also be linearized, via separate 
control 4a¥> ^loops according to the 'modular' approach to linearization disclosed as part of this 
invention. 

All of the effects mentioned above vary with time and circumstances. They are nonlinear 
and thus distort the sound wave shape, in both amplitude and phase, relative to the input audio 
information. Furthermore, due to the inherently bi-directional nature of the transconductance and 
the electromechanical transduction, and of the coupling between them, distortions in any one 
process ar e mirror e d in can affect any of the other processes. Most importantly, it is the 
nonlinearities inherent in the electromechanical transduction whieh- that make the linearization 
and control of the overall process very difficult in prior art. 

Page 23, line 8, please amend the following paragraph: 

There are well-recognized nonlinearities in the drive current as a^function of voltage, 
caused by the dependence of effective coil impedance and of the motor's BEMF upon coil 
position relative to the magnet assembly. The effective spring stiffness of the coil/diaphragm 
assembly, likewise dependent on coil position, as is the motor factor, result in well-recognized 
sources of nonlinearity. Additionally, more gradual changes of coil impedance due to Ohmic and 
environmental heating cause the drive-current response to vary over time. All these effects cause 
power and frequency dependent distortions of the audio signal. 

Page 24, line 9, please amend the following paragraph: 

The three processes can be described by a mathematical model, comprising a system of 
coupled equations specifying the rate of change (evolution) of each of a complete set of state 
variables, such as coil current and coil position, at any given time, in terms of the state vector at 
the same and all previous times. Such equations are termed "integro-differential equations", and 
are nonlinear in the case at hand. In the prior art, the model equations are usually approximated 
as having no "memory", in the sense that the rates of change of state variables are taken to be 
wholly determined by (generally nonlinear functions of) state variables at the same instant of 
time; such memory-less evolution equations are simply termed "differential equations". The 
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mathematical model according to th e pr e sent invention, however, includes m e mory e ffects, as it 
has been discovered that they cannot, in general, be e ntirely neglected in mod e ling th e audio 
reproduction system. 

Page 24, line 32, please amend the following paragraph: 

A nonlinear process can be very complex, and the number of terms kept in the evolution 
equations, as well as the decision whether or not to include memory effects, and if so which 
ones, can vary depending on the degree of approximation required in the control methodology. In 
the explanation whieh -that follows, it will be seen that simplifying the approximations to the 
most basic mechanisms of the three processes yields several coupled "ordinary" nonlinear 
differential equations. Anyone skilled in the art will appreciate that using approximations is a 
compromise, and that beyond a certain point, enlarging or truncating the list of modeled effects 
does not alter the fundamentals of the invention. 

Page 27, line 21, please amend the following paragraphs: 

Most of the parameters and parameterized functions appearing in equations (6) through 
(4410), specifically R^, R^, £l(x), K{x), h{t) and the functions gj through depend on 
temperature, which is assumed to vary slowly as compared with timescales characterizing audio 
response. For the approximation to be fiilly self-consistent, the acoustic-load part of R^ should 

actually be replaced with a memory term related to h{t) ; the fact that a constant R^^ is instead 
used in equation (9) is a further, non-essential approximation. 

The time integrals in equations n\ and (8) encode memory effects due to eddy currents, 
while the integral in the pressure equation (10) encodes memory effects due to acoustic 
reflections and dispersion. All of these integrals represent the dependence of the rate of change 
of state variables at any given time, upon the history (past values) of those same state variables. 
Although effects from the infinitely remote past are in principle included in these integrals, in 
practice the memory of past positions and currents fades eventually, because the audio signal is 
band limited. 

It has been found that, while the memory effects encoded in equations (7), (8 Vand (10) 
are important for modeling the dynamics of an audio reproduction system, they are second ■ 
efder of secondary importance in the context of a distortion d i stort ion- correction controller. 

The spectral contributions to the dynamic coil excursion x{t) are dominated by low 
frequencies, a fact well recognized in prior art. In consequence, it is often a reasonable 
approximation to replace the delayed positions x(r), x(r,) and x{t^) in the memory integrals 

-10- 
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of equations (7)-(8) with low-order Taylor expansions about the present time (i.e. about r = / , 
r, = t and, = / respectivelyT)^ In this way, positional memory effects are neglected^f while 
the more important memory effects involving delayed response to current and velocity, are still 
included. If this further approximation is implemented, and terms quadratic and higher in coil 
velocity are neglected, the electromechanical and elastic parts of the above system of evolution 
equations, equations (6) through (9), simplify to the following form. 

The coil-circuit equation (governing the transconductance component of the signal 
conditioning process) becomes: 

v^At) = Re m + ^(00^^,(0 + r^,,,(0 (1 1) 

— Wh^where now <I>^_^„^„.,(0 and F,^,/^(0 simplify to 

^Jyna^^c (0 = Bl(x(t)) + jdv g, - T,X(t)) /(Z") + 

(12) 

IT, ^ ^ 

+ fdz, \dT,g,^''\t-T„t-T„x(t))i(T,)i(T,) 

—00 —00 

— and 

/ / r, 

—00 —00 — 00 

f 

+ jdTg,it-T,xit))xiT)iiT) 

-00 

(13) 

— respectively. 

Page 30, line 4, please amend the following paragraphs: 

— Whefe -where is a constant. Since all memory and eddy-current effects have been 
suppressed in equations (14)-(16), parameter estimation of L^(x)^R^ and from empirical data 
will show that they are frequency-range dependent; and, that, furthermore, actually depends 
upon x(t) since it includes the resistive counterpart to effective coil reactance L^{x) caused by 
eddy currents. 

Equation (14) is an oversimplification. As recognized in the audio industry, a transducer 
voice coil is characterized by a frequency-dependent complex effective impedance, which we 
denote (o) , x) to indicate that it also depends upon coil position; it also implicitly depends 
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upon other, more slowly varying parameters, such as temperature. The effective coil impedance 
Z^{q),x) characterizes one aspect of the relation between voltage signal V^^.f{t) applied to the 
voice-coil circuit on the one hand, and the coil current /(/) caused by this voltage, on the other. 
This voltage-current relation, or functional, as it is known mathematically, is nonlinear, and 
furthermore involves electrodynamical memory effects (distributed delays) as described above. 
In general this relation can be expanded in a functional series of the type known in the literature 
as a Volterra series. The multivariate coefficient-functions of this Volterra series depend on coil 
position and motion within the magnetic-circuit airgap. 

Current-nonlinear effects, i.e. deviations from linearity of the voltage-current functional, 
were found to be measureable. For the Labtec Spin70 speaker transducer, one of the large signal 
data parameters which ar e illustrated in FIG. 4, namely , was found to vary with i{t) as the 

coil neared its negative excursion. However, it was also found through modeling, simulation and 
measurements that current-nonlinear effects in speakers are typically small (at the few percent 
level), although they can become important for woofers played at high volumes. Thus, for many 
transducers, the fiill complexity of the current response /(/) to a given applied voltage V^^^i{t) 

can often be usefully approximated by a linear functional relation, in which memory effects (due 
to eddy currents in the magnetic pole structures, and in the aluminum coil former if any) are still 
included. This approximate linear relation can be derived from equations (1 1)-(13) and is 
expressed as follows: 

V,^{t)'=^RJ{t)\v{t)Bl{xit))\ \dT g,{t T,x(t))i(T) 

-00 

/ 

v^,{t)=R,m+mBm))+ g,{t-T,x{tmT) (i?) 

-00 

Page 31, line 20, please amend the following paragraph: 

The second (velocity dependent) term on the right hand side of equation (17) is the 
BEMF due to coil motion; the other two terms comprise the EMF due to the overall effective coil 
impedance. Within the approximation, invoked above, of a slowly changing (low frequency) 
position x{t) , the Fourier transform of with respect to time is simply the subtracted effective 
coil impedance in frequency domain, i.e. the coil impedance with the Ohmic coil term 
subtracted. We denote this subtracted coil impedance as Z'"^e{(o,x) . More precisely, when a 
probe voltage signal at a typical audio (or supersonic) frequency is applied to the voice coil and 
the attached diaphragm is mechanically held (blocked) at a fixed position x , the effective 

-12- 
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impedance, due to the coil's inductance and its interaction with eddy currents and magnetization 
within the magnetic poles, is by definition- Z^{co,x) = Z'''\{o),x) + , where the term is 
added in series and represents the coil's Ohmic resistance (see Equation (17)). Note that the 
subtracted impedance Z^"*c(<a;, jc) has both resistive and reactive components; the former is 
attributable to eddy-current dissipation inside the magnetic poles (and also in the coil former, in 
case that is made of aluminum). The reactive component of Z"'*e(fi;, x) is known in prior art as 
L^ix), with the frequency dependence often left implicit, as it was in equations (14)-(15) above. 

Page 32, line 18, please amend the following paragraph: 

For sufficiently high frequencies, and in the case of non-metallic former, the subtracted 

sub 

impedance (co, x) arises from currents and EMF's induced in the coil and within a narrow 
skin layers, within the pole structures and adjacent to the coil. For a simple cylindrical geometry 
with infinite axial extent, Z'"^e(co,x) is independent of jc ; in that approximation, Vanderkooy 
[J. Vanderkooy, J. Audio Eng. Soc, Vol. 37, March 1989, pp.1 19-128] -has shown that the 
(complex plane) phase angle of the subtracted impedance begins to approach an asymptotic 
value of 45'' once the frequency increases well above the normal modes of mechanical 
resonances. Measurements for actual speaker transducers yield a range of possible asymptotic 
phase angles, both above and below this value [J. D'Appolito: "Testing Loudspeakers", Audio 
Amateur Press; 1998.] For the tabTeeLabtec Spin 70 speaker transducer analyzed in the present 
study, the asymptotic phase angle was measured to be approximately 70" , varying little with 
coil/diaphragm position x . 

Page 33, line 28, please amend the following paragraphs: 

(II) ^ — Progressively smaller nonlinear effects can be corrected by applying successive 
new linearizing filters, and this progression of successive corrections will often converge in the 
sense of perturbation theory. 

It should be noted that the ability to systematically apply more and more modular control 
tiers^ can be useful even if a higher-tier correction is larger than a lower-tier one. 

FIG. 7 is a flow chart whj ^that illustrates the process of linearization in accordance with 
the present invention. First, a model of a portion of the audio reproduction system is provided in 
step 1301. Next, a control engineering technique is utilized in the time domain to control an 
output of the audio transducer system based upon the model, via step 1302. 
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Page 34, line 20, please amend the following paragraphs: 

FIG. 8 is a block diagram of the main portion of a sound reproduction system and a 
control system for controlling the operation of the sound reproduction system in accordance with 
the present invention. An audio signal 1401 is input to a controller 1402, which contains 
algorithms based on a control model, which in turn is based on a physical model (such as the one 
described by equations (6)-(16) of this section) of the processes within the audio transducer 
system. These algorithms may be functions of state variables such as acceleration, velocity, and 
position of the coil/diaphragm assembly. With reference to FIG. 6, the modeled processes may 
include the signal conditioning process 1 102, the voice coil transduction process 1 104, and the 
sound conditioning process 1 106, as discussed above. The state variables 1403 from the sound 
reproduction processes are input to the controller 1402 from a measurement system 1404. The 
measurement system 1404 consists of a sensor conditioner 1405 and a plurality o f one or more 
sensors, 1406a, 1406b, and 1406c, which take measurements of variables from the sound 
reproduction system. The sensor conditioner 1405 amplifies and converts the signals from the 
sensors 1406a, 1406b, and 1406c to the state variables 1403, which are provided to the controller 
1402. Sensor 1406a may, for example, measure a variable such as current from the drive 
amplifier 1407. Sensor 1406b may, for example, measure an internal circuit parameter, such as 
parasitic capacitance, of the transducer 1408. Alternatively, sensor 1406b could electronically 
measure the impedance of one of the voice coils of transducer 1408, or it could optically measure 
an indicator of voice coil position. Sensor 1406c may, for example, measure a variable from the 
► acoustic environment, such as sound pressure by using a microphone. By digitizing both the state 

variables 1403 and the audio signal 1401, and combining them via a DSP, the controller 1402 
modifies the audio input 1401, converts it back to an analog voltage, and thus outputs a 
compensated analog audio signal on line 1409 to the amplifier 1407. The amplifier 1407 outputs 
a drive signal on line 1410 to the transducer 1408. 

The audio transducer state variables w^neb -that are measured and fed back to controller 
1402 are generalized coordinates of the transducer dynamical system. These generalized 
coordinates usually vary nonlinearly with the position of the voice coil/diaphragm assembly with 
respect to the transducer frame, and thus, with suitable calibrations, serve to provide controller 
1402 with estimates of recent values of that position. Controller 1402 then uses these real-time 
position estimates to suitably modify the input audio voltage signal before applying it across the 
voice coil. Multiple position-indicating signals can be fed to the controller, as depicted in FIG. 
8; they are derived from one or more position-indicating generalized coordinates. It may be 
useful to measure more than one position-indicating generalized coordinates, because in some 
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portions of the range of coil/diaphragm excursions, it could happen that a given generalized 
coordinate may not be a monotonic function of coil/diaphragm position, while another 
generalized coordinate is monotonic in that portion of the range. Thus, the advantage of 
measuring and feeding back values for multiple generalized coordinates, is that these coordinates 
may be chosen in such a way that the configuration space of their joint values is approximately a 
one dimensional differentiable manifold, where the coil/diaphragm position is a continuous and 
differentiable function on this manifold. And if each of the selected generalized coordinates is 
also a continuous and differentiable function of coil/diaphragm position, the mapping between a 
tuple of simultaneously measured generalized coordinates and the corresponding position, is both 
invertible and differentiable, allowing the use of the tuple to compute the audio signal 
modification within the controller DSP. One embodiment of this computation, based on a single 
generalized coordinate whieh -that is derived from infrared optical measurements, is described in 
detail in Detail 10. 

It will be readily apparent to those skilled in the art that additional and different sensors 
may be utilized, and different signal conditioners may be used to recover state variables and 
internal parameters from the sensor signals and provide control signals to the system. Additional 
sensors may include, for example: accelerometers, additional transducer coils, or new coil-circuit 
elements. Such sensors can provide analog measurements of various voltages appearing in the 
transconductance equation (14), or of other voltages whieh -that allow the estimation of various 
terms and state variables in either equation (14) or the mechanical (transduction process) 
equation (15). State variables and parameters must be identified for each of the sound 
reproduction processes, and a sufficient set of them must be measured to effect control. 

It has been discovered that measurements not usually regarded as state variables can be 
used effectively in controlling the audio reproduction processes. In the prior art systems, the 
following variables are typically considered as defining state: 

X axial position of coil/diaphragm assembly, 

X axial velocity of coil/diaphragm assembly, 

jc axial acceleration of coil/diaphragm assembly^ 

/ voice-coil current. 

What follows is a list of other measurable variables, among them internal parameters 
characterizing the processes whieh -that are considered constants in small signal analysis, as well 
as state variables, such as pressure, which would be externally measured (using a microphone in 
this case). The variables and parameters on this list can all be used in practicing the present 
invention. Control systems using one or more of these variables and parameters are described 
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below. Some measurable variables can be measured by reference to other variables through 
known functional dependencies; for instance, temperature can be inferred from coil resistance 
and a lookup table. Internal -parameters and other variables not listed in the above list include, 
for example: 

V(t) voice coil voltage, 
/(/) voice coil current, 
voice coil resistance, 
voice coil inductance, 
complex voice coil impedance, 
C parasitic ^oicc coil/magnct parasitic capacitance, 
BEMF back-EMF, 

^ complex phase angle of voice coil impedance, 
voice coil temperature. 

There are other internal parameters such as Bl and K, respectively the motor factor and 
suspension stif&iess. These parameters may be difficult to measure directly, although they can be 
extracted from measurements of other variables via parameter estimation methods. The voice- 
coil voltage V{t) and voice coil current^ /(/) are considered as-internal variables, rather than 
stimuli, because the full audio transduction process according to the present invention includes 
creating V(() and /(/) as internal variables. 

Page 39, line 13, please amend the following paragraphs: 

The control model treats the motor factor Bl{x) , the effective coil inductance L^(x), and 
the suspension stiffriess K(x) as functions of jc(/) , the current axial position of the 
coil/diaphragm assembly. These three functions cause most of the nonlinearities, and thus 
distortions, of audio transducers, as explained above. The motor factor Bl(x) determines the 

motive force term in equation (19) as well as the BEMF term in equation (18); L^{x) determines 
the inductive EMF term in equation (18); while K(x) determines the elasto-acoustic restoring 
force in equation_(19). In the context of the present invention, these three functions are derived 
from calibration measurements on the system, which yield the functional dependence of Bl , 

and K upon x ; these functions can, for instance, be obtained from commercially available 
transducer test equipment such as a Klippel GMBH laser metrology system. In one embodiment 
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of this invention, the functional dependences Bl(x) and L^{x) are entirely obtained from such a 
laser metrology system, while K(x) is obtained by combining knowledge of Bl(x) and L^ix) 
with ramped DC-drive calibration runs, as fully described in Details 5 and 10 below. 

In transducer operation, the three functions Bl(x) , L^(x) and K{x) must be combined 

with approximants to a function mapping the measured position-indicator state variable onto the 
actual position x, as described in Details 4,.5, and 10 below, in order to provide the controller 
DSP with an estimate for the values of Bl , and K{x) at the current present moment t . 

The controller then estimates the BEMF term by multiplying the estimated present value 
of Bl{x{t)) by an estimate for the present velocity x{t) ; the latter may be obtained either from a 
numerical differentiation of the recent history of discrete position measurements, or from an 
independent velocity measurement. In one embodiment of the present invention, velocity is 
estimated via numerical differentiation of estimated position, as described in Detail 10 below. 
Simulations of the BEMF correction shows that it can be usefully filtered in the frequency 
domain, as this correction has its greatest effect over a limited frequency range. Such filtering 
reduces the noise due to the numerical differentiation of position. Once the nonlinear BEMF 
term Bl{x) x in equation (18) is thus estimated, it is corrected for by being added by the control 
circuit to the voltage representing the audio information. A linear BEMF term can also be 
calculated and subtracted from the voltage representing the audio information, in order to 
provide damping if required. The subtracted linear part of the BEMF is chosen such that the 
effect of the subtraction is to electronically add back a positive constant to the mechanical drag 
coefficient in equation (19). This positive constant is some adjustable fraction, p , of the 

Thiele-Small small-signal BEMF contribution to the drag coefficient that would arise due to the 
equilibrium value 5/(0) m^/7/iow/ any correction. Thus, 

In many cases of interest the effective coil inductance {x) in equation (18) is very 

small. If we neglect this inductance, the inductive EMF term L^(x) — in equation (18) 

dt 

disappears, and that differential equation becomes an algebraic equation. With this 
simplification, the voltage signal that is output from the control circuit to the voice coil in order 
to compensate for the nonlinear BEMF is: 



_ 5/(0) J . , Bl{0) 



(20) 
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A\^hefe- where V^^^^^.^ is the voltage representing the audio signal before the BEMF correction. 
Note that other modular corrections may be included in V^^j^^ , as described below. 

We next turn to the case where the effective coil inductance LJx) in equation (18) is not 

neglected, and describing describe another type of modular control law in the context of the 
present invention, namely a control law that correcting corrects for the inductive EMF term in 
equation (18). Like the BEMF control law described above, the inductive control law partially 
linearizes the transductance sub-process. Specifically, the inductive control law addresses the 
nonlinearity, and thus distortion, caused by the position dependence of the effective coil 
inductance L^{x) . In order to derive the inductive control law in as simple a manner as possible, 

the BEMF term is temporarily ignored in the transconductance equation (18); later in this 
section, all four of the modular control laws described in the context of this invention (BEMF, 
inductive, spring and motor factor) will be combined. 

Since the embodiment described below for the correction of the inductive EMF term 

L^(x) — in equation (18) has no history in prior art, the derivation of this correction is 

presented in some detail here. For simplicity, noise is ignored in this derivation, as is-are the 
deviations of the in-operation digital signal processor (DSP) estimates for L^{x{t)) from the 
actual values of this variable. 

Page 42, line 5, please amend the following paragraphs: 

If f^a„rf/o(0 9 ^(0 treated as known functions, equation (23) can be viewed 

as a linear first-order ordinary differential equation for the unknown function i(t) . It is a welU 
known mathematical fact that this differential equation admits a unique solution i{t) for any 
given causal signal V^udioiO ? i-^- for an audio input signal that begins at some definit e initial time 
tQ in the pas t given an initial condition i(t^) = Iq . Th e latt e r condition can saf e ly b e assum e d, 

sSince any real-life signal is causa K we can safely assume that there is an initial time such that 
/(^^) = 0 and V^^^iJip) = 0 . On th e oth e r hand Then, it is easily verified by substitution that a 
particular solution of the differential equation (23) is given by 

m = V.:,MfK (24) 
The combination of these two facts, namely, that equation (23) has a unique solution for the coil 
current in terms of the audio voltage input, and that equation (24) is a particular solution of 
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• 

equation (23) and is valid at an initial time such that Ht^) = 0_and_F^^^.^(/o) = 0^— completes 

the proof that equation (24) does in fact hold for all values of / . In other words, it has been 

proven that the coil current i(t) is related to the audio signal V^^^^t) by a simple Ohm's law, 

without any inductive term, provided that BEMF is ignored and that the control law of equation 
(22) is implemented. 

This demonstrates that by simply adding to the audio signal voltage a term that is the 
derivative of this same audio signal, multiplied by the ratio of the nonlinear inductance to the 
coil resistance, as done in equation (22), a correction for the effects of inductance alone can be 
made. In one embodiment of the present invention, the voltage differentiation on the right hand 
side of equation (22) is implemented numerically by the DSP, as fully described in Detail 10 
below; this alone introduces additional terms on the right hand side of equation (24), thus 
making the elimination of the inductive term approximate, rather than exact-. Furthermore, it will 
be appreciated from the detailed description of polynomial interpolations in the context of this 
invention (Detail 10 below) that the correction of the inductive effect by the physical controller, 
as opposed to the ideal one assumed in the above derivation, is approximate, rather than exact. 
This caveat would hold even were an exact, analog differentiation to be used by the controller. 
And it also holds for the numerical BEMF correction described above. 

In the case of input to a voice coil which that is used for audio reproduction, removing all 
the inductance as described in equations (21)-(24) might lead to an equalization problem, since 
the higher frequencies can be over-compensated. Thus, in one embodiment, an optional linear 
part of the inductance is added back to endow the audio system with a flatter frequency response. 
This is described in Detail 10 below. 

Page 43, line 32, please amend the following paragraphs: 

The correction of the nonlinear electromechanical effects in the mechanical 
(transduction) equation of motion (19) is based upon a derivation similar to, but different from, 
the standard control theory derivation of a control equation presented in the Background section 
above as prior art. One practical problem with the mechanical equation (19) as a starting point 
for a control models is that the inertia term involves the coil/diaphragm acceleration jc . This term 
increases rapidly with frequency, eventually becoming too large to be considered in a 
compensation system. However, because the acoustical radiation efficiency of the cone also 
increases with frequency, the inertia non-compensation is balanced by the radiating efficiency, 
within limits. This trade-off is known in prior art to result in a more or less constant output over 
a range of fi-equencies referred to as the 'mass controlled' range. Transducers are normally 
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designed with this effect in mind. 

By ignoring mass in equation (19), that is to say by neglecting inertial effects, the 
following first order differential equation is obtained: 

K,x^xK{x)=Bl{x)i{t) (25) 
In the general nonlinear state space form, equation (25) is recast thus: 

X = ^{x) \ if/{x) u{t) x= (p^ (x) + ij/^ {x) u{t) (26) 

— where, 

./ ^ -xK{x) ^ -xK{x) /X Bl{x) /X Bl(x) 
^ms ^ms ^ms ^ms 



(27) 



— and: 



u{t)=i{t) (28) 
Following the feedback linearization approach, consecutive derivatives of the transducer 
output are taken until its input, u(t) , appears in one of the derivatives. But that is already the 
case in equation (26), which^ when combined with equation (27), yields for the first derivative of 
coil/diaphragm position x{t) : 

,^ -xK{xhBl{x)u{t) ^29) 

Note that the input, u{t) , indeed appears explicitly in the first derivative of the position 
state variable, x . 

The controller linearizing the transduction process should cause the transducer output 
x(t) to be proportional to the audio input. Equating x{t) with V^^^.^{t) in equation (26) and 

solving for w(r), and assuming that the function tff{pc)-i//^{x) defined in (27) is nonsingular, we 

obtain: 

u{t)A^{^))r[-M ^ r4{t) = [if^, (x))r [- {xh w] 

(30) 

— ^where w(t) is the generator or reference (in our case the audio program input V^^^^^{t) to the 

uncorrected transducer), and u(t) is the actual voltage input to the voice coil in the controlled 

(corrected) transducer if the signal conditioning process is ignored. Substituting and rearranging 
terms in equations (27),_(28) and (30), provides: 
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By applying this (ideal) control equation to the second order differential transduction 
equation (19), it is possible to see whether the latter is thereby linearized. 
Substituting equation (31) into equation (19) provides: 



mx I Rms^ I K(x)x = Bl{x) 



Bl{x) Bl{x) 



mx + R^x + K(x)x = Bl{x) 



xK{x) 

Bl{x) ' " Bl{x) 



+ w- 



(32) 



This leaves, 

^^ + Ks^ = ^Rms (33) 

Equation (33) is a linear differential equation with constant coefficients. Note that from 
the above a general method of linearizing this form of nonlinear dynamical equation is presented, 
and any further linear terms can be added to the equation without changing the validity of the 
linearization approach. 

Lumping the terms of the rearranged control equation (31) and using equation (28) 
provides the following form of the transduction control equation: 

u{t)^S{x) + w{t)B{x) (34) 
— Whefe -where S(x^ and B{x) are functions of position and w{t) is the audio information. 

Equation (34) provides a correction for the open loop non-linear transfer function of the 
speaker transducer, provided that the dependencies of S{x) and B{x) on x are known and that 
real-time measurements or estimates of x are made available to the controller during transducer 
operation. 



Page 46, line 14, please amend the following paragraph: 

Clearly, the control law given by equation (34) removes all restoring force due to the spring; 
a thus corrected transducer would not be stable. Thus a linear (non-distorting) restoring force 
must be subtracted from jc K{x) . The magnitude of the effective spring constant of this residual 
electronic linear restoring force^ can be selected based on the required resonant frequency. This 
then in effect reduces the transducer operation to the linear case of zero motor factor and a linear 
(Hooke's law) elastic restoring force. A full description as to how this subtraction is 
implemented in one embodiment of the present invention, is presented in Details 5 and 10 below. 

Page 46, line 32, please amend the following paragraphs: 

The control model of equation (34) applies only to the transduction process itself; i.e. it is 
based on a model of the current to velocity transduction process, and does not cover the process 
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of injection of current into the coil (the signal conditioning process); nor does it cover the 
radiation of the sound waves out of the speaker enclosure into the acoustic environment (the 
sound conditioning process). Likewise, the control models of equations (20) and (22) above, 
suitably combined, eliminate or reduce only those nonlinearities arising from the 
transconductance component of the signal conditioning process, but do not correct either of the 
other two processes (transduction or sound conditioning). And all of the above control laws can, 
and have been, applied together, or in various partial combinations, in the context of the present 
invention. This illustrates the modularity of the control approach described as part of the present 
invention, as discussed in Detail 1 above. Furthermore, the transduction control law of equation 
(34) can be subdivided into "spring correction" and "motor factor" modular units; e.g. if only the 
first term on the right hand side of equation (34) is used, this represents a control law which only 
linearizes the elastic restoring force. Thus, the number of modular control laws described by the 
above equations can actually be counted as four: BEMF, inductive, spring, and motor factor. 

If a choice is made to simultaneously implement all of these modular corrections: the 
BEMF correction (equation (20)), the inductive correction (equation (22)), and the transduction 
corrections (equation (34)), this can for example be done as follows. The last term B EMF 
corr e ction of equation (20) is added to the voltage given by the right-hand side of equation (34); 
and-then the new overall voltage, u^(t) , still in the digital domain, is numerically differentiated 
(as described in Detail 10 below), and this numerical derivative is fewHy-combined with 
w,(r) itself in accordance with equation (22). Finally, the BEMF correction term of equation (2) 
is added to the new voltage. The overall combined control model for the coil voltage is thus as 
follows: 



u{t) = u^ (0 I 




(36) 



where 



on 



As explained above, the precise order in whicli the modular corrections are applied is not 
very important, as has in fact been demonstrated in the context of this invention. 
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In order to add back an effective electronic linear restoring force, as discussed above and 
in Detail 5, the term S(x) on the right-hand side of equation (35) must be replaced by the 
subtracted version, 

^ ^ Bl(x) 

where q is the fraction of the uncorrected suspension stiffness at equilibrium that is added back 
electronically. Thus equation (35) now becomes, 

«,(0 - S{x) - ^^^^Q) I wBix) + iBlix) - pBliOy/Blix))x(t) ^ 
Bl(x) 

^u,it) = S{x)-^^ + wBix)- pR^^m , (39) 

Bl (x) Bl{x) 

while equation (36) remains unchanged. 

In case a choice is made to implement only the transduction correction law, it is still 
necessary to perform the suspension stiffness subtraction, for stability purposes - as explained 
above. Thus, the flill transduction control law in accordance with the present invention is the 
following modified version of equation (34): 

u{t) = S(x) - + w B{x) (^40) 

Bi (X) 

One view of the control method described in this invention is that it belongs to the genre 
of feedback linearization controllers. The transconductance component of the signal conditioning 
process, and the transduction process, together may be thought of as a dynamic system with 
voltage input and displacement output. The dynamics of this system are governed by a physical 
model that can be represented as a three-state system with current, displacement, and velocity as 
its state variables. As seen above, despite the interactions among all processes comprising the 
audio reproduction system, various processes and sub-processes can be separately controlled 
according to this invention by applying only one of the separate basic linearization control laws 
encoded by equations (20),(22), and (34), or these control laws may be applied in various 
combinations - depending on user preferences. One option is to apply all of them, as encoded in 
equations (36) and (3739), as well as in equations (56) (59) (6n-(64) in Detail 10 below. 

Page 49, line 14, please amend the following paragraphs: 

FIG. 9 shows the feedback linearization process 20400 with the control law of equation 
(34), which only linearizes the transduction component of the signal conditioning process, 
without an electronically restored linear restoring force. The audio signal, V^^^.^ = w 20401, is 
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input to a Linear Compensation Process module 20402 (henceforth abbreviated as LCP). The 
LCP 20402 multiplies w by the compensation function B(z) , where z 2041 1 is the estimated 
present value of the position variable. The present value of position variable z 2041 1 is obtained 
from the transduction module 20408 of the three-state overall transducer system, via a two step 
process: first the position indicator state variable f{x) 20413 is measured by the positional 
sensor module 20412, and then the value of f(x) 20413 is fed as input to a sensor inversion 
module 20414, which estimates actual position x via an interpolation method as described in 
Details 5 and 10. Actual position x 20409 and actual velocity x 20410 are fed from the output 
of transduction module 20408 back into the input of the transconductance module 20406, via the 
physical system itself (not as measured data). The estimated x value, z 204 11, is inputt e d fed 
into to the LCP 20402 and also to an S -lookup module 20415. The output of module 20415, 
5(z) « Six) 20416, as well as the LCP output B{z)w 20403, are both fed as inputs to a 
s umme r summing junction 20404, the output 20405 of which is the corrected audio signal ( V^^-j 
of equation (34)). This corrected audio signal 20405 is provided as input to the transconductance 
module 20406 ofthe three-state transducer system. The current output 7^^^, 20407 of the 

transconductance module 20406 is provided as input to the transduction module 20408. 

FIG. 10 shows the feedback linearization process 20500 for the control law given by 
equation (3^40); again only transduction corrections are made, but now a linear spring constant 
(suspension stiffness) is electronically added, as explained above and in Detail 5. The audio 
signal, V^^^.^ = w- 20501, is input to an LCP module 20502. The LCP 20502 multiplies w by the 
compensation function B(z) , where z 20514 is the estimated current value of the position 
variable. Value z 20514 is obtained from the transduction module 20508 of the three-state 
overall transducer system, via a two step process as in FIG. 3049: the positional sensor module 
2051 1 outputs the measured position indicator state variable f(x) 20512, and measured state 
variable /(jc) 20512 is fed as input to a sensor inversion module 20513, which estimates actual 
position X via the interpolation method. Actual position x 20510 and velocity x 20509 are fed 
back from the output of the transduction module 20508 to the input of the transconductance 
module 20506 via the physical system itself 

The estimated x value, z 20514, is this time inputt e d to fed into three modules: to the 
LCP 20502, to an S -lookup module 205 16, and to a new 'Electronically Restored Linear 
Spring'(henceforth ERLS) module 20517. The output of module 20516, S(z) ^ S(x) 20415, as 
well as the LCP output B(z)w 20503 and the output 20518 ofthe ERLS 20517, are all fed as 
inputs to a summe r summing junction 20504, the output 20505 of which is the corrected audio 
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signal (V^^.f of equation (34)). The corrected audio signal 20505 is provided as input to the 

transconductance module 20506 of the three-state transducer system via the physical system. 

FIG. 1 1 shows the feedback linearization process 20600 for the control law given by 
equation (^39) alone, without the inductive correction (36); i.e. for a control law correcting for 
spring, motor factor and BEMF nonlinearities, including an electronically restored linear spring 
and electronically restored contribution to the linear drag force term, as explained above. The 
audio signal, = w 20601, is input to an LCP module 20602, The LCP 20602 multiplies w 
by the compensation function B{2) , where z 20622 is the estimated present value of the 
position variable. The output B(z)w 20603 of the LCP module 20602 is provided as input to the 
summe r summing junction 20604. Value z 20622 is obtained from the transduction module 
20610 of the three-state overall transducer system, via a two step process as in the previous 
figures: the positional sensor module 20613 outputs the measured position indicator state 
variable f(x) 20614, which is then fed as input to a sensor inversion module 20615. Sensor 
inversion module 20615 estimates actual position jc via the interpolation method. And as in 
previous figures, the actual position x 20612 and velocity jc 2061 1 are fed back by the actual 
physical system from the output of the transduction module 20610 to the input of the 
transduction transconductance module 20608. The estimated x value, z 20622, is now inputt e d 
te fed into four modules: to the LCP 20602; to the S -lookup module 2061 8; to an ERLS 
module 20620; and finally, to a BEMF-computation module 20616, which applies a numerical 
differentiation operation D to z 20622. The output 20619 of the module 20618, as well as 
output 20621 of module 20620 and output 20603 of the LCP 20602, are summed in the 
summ e r summing junction 20604. The output 20605 of summe r summing junction 20604, along 
with the output 20617 of the BEMF-computation module 20616, are provided as inputs to a 
second summ e r summing junction 20606; finally, the output 20607 of the second 
summe r summing junction 20606 is the corrected V^^^, , which is provided as analog input to the 
transconductance module 20608 of the three-state transducer system. And the analog coil current 
I^^^, 20609, output by the transconductance module 20608, is provided by the physical 

transducer as input to the transduction module 20610. 

FIG. 12 shows the feedback linearization process 20900 for the control law given by 
equations (36) and (3^39), i.e. implementing all the corrections described in this section, and also 
implemenfing two numerical Low Pass Filters: one between the position-indicator variable 
measurement and the sensor inversion, and another after the computation of the fiilly corrected 
coil voltage and before it is fed as input to the coil. The audio signal, V^^^.^ = w 20901, is input to 
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an LCP module 20902. The LCP module 20902 multiplies w by the compensation function 
B{Zj^), where Zj- 20921 is a filtered version of the estimated present value of the position 

variable. The output -B{2)w- B{z y^)w 20903 of the LCP module 20902 is provided as input to the 

summe r summing junction 20604. Value Zj^ 20921 is obtained from the transduction module 

20910 of the three-state overall transducer system, via a three step process: the positional sensor 
module 20912 outputs the measured position indicator state variable f{x) 20913, which is then 
fed as input to the low pass filter LPF2 20924, the role of which is to suppress sensor noise; 
LPF2 would typically roll off at 1-2 kHz. The output 20925 of LPF2 20924 is fed to the sensor 
inversion module 20914. Sensor inversion module 20914 again estimates actual position jc via 
the interpolation method, in the digital domain; while the actual position x 2091 1 and velocity 
X 20912 are fed via the physical transducer plant, back from the transduction module 20910 to 
the transconductance module 20908. The estimated jc value, now called Zj. 20921,- is inputt e d 

t efed into the following three modules: to the LCP 20902, to the ERLS module 20920, and to the 
BEMF-computation module 20915, The S -lookup module 20917 receives its input this time 
from the filtered, but not inverted, positional indicator variable measurement result 20925. The 
outputs of the four modules 20915, 20917, 20919 and the LCP 20902, labelled respectively 
20916, 20918, 20920 and 20903, are summed in the Sttmrne rsumming junction 20904. The 
output 20905 of summ e r summing junction 20904 is passed to an inductive-correction module 
20927, which again applies a numerical differentiation operation D , this time to the numerical 
output voltage 20926 of the summ e r summing junction 20904. The output 20906 of the 
inductive-correction module 20927 is provided along with numerical output voltage 20926 
multiplied by to a second s umm e r summing junction 20928, whose output- 20907 is fed to the 

low pass filter LPFl 20922. The low pass filter LPFl 20922 implements a (partial) correction for 
the voice coil inductance at equilibrium. The output 20923 of LPFl 20922 is finally fed as the 
corrected analog voltage F^^^, to the transconductance module 20908 of the three-state transducer 
system. As in the previous figures, the physical transducer plant provides the analog output 
current /^^./ 20909, output by the transconductance module 20908, as input to the transduction 
module 20910. 

As emphasized above, the present invention requires at least one state variable to be 
measured in operation for any given run. In the control diagrams depicted in FIG. 9, FIG. 10, 
FIG. 1 1 and FIG. 12, it has been assumed for convenience that only a single state variable is 
measured (although at least two variables, such as for example the position measured by an IR 
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sensor x-^ and the position measured by a laser sensor x^.^ ^ x , would need to be nneasured 
during offline calibration runs in order to derive an interpolated function / (x)). 

The process of applying a state variable feedback law based on a plurality of 
measurements of one or many state variables^ is depicted in FIG. 13. The process 21000 begins 
with one or several measurements of a state variable or variables from a plurality of sensors, 
21001 through 21002. For example a transducer's coil/diaphragm displacement, x , may be 
measured both via the parasitic capacitance method (Details 7 and 12 below) and the IR method 
(Details 8 and 13 below). The respective state variable measurement signals, 21003 through 
21004, are passed from the sensors to the state estimation module 21005, which synthesizes the 
desired partial or full state variable estimate, 21006, which in general is a vector state variable. 
This state variable estimate 21006 is in turn used in the application of the control law 21007 in 
place of the actual state variable. 

For all practical purposes, none of the sensors, 21001 through 21002, can measure its 
intended state variable exactly. The measurement is always corrupted to some extent by factors 
including nonlinearities in the measurement, measurement noise, quantization noise, systematic 
errors, etc. The task of the state estimation module 21005 is to mitigate these corrupting effects. 
This task may include all or some of the following ingredients: inverting the nonlinearities of the 
sensors to provide a more linear response to the measurements 21001 through 21002; adaptation 
to minimize the sensitivity of the state variable estimate 21006 to parametric uncertainties in the 
measurement, such as uncertainty in gain; filtering the measurement signals 21003 through 
21004 to minimize the effects of noise; or fusing multiple measurements of a state variable into 
one state variable estimate 21006. In addition, many engineering objectives are taken into 
consideration in the design of the state estimation module 21005. The tradeoffs include such 
d e sid e rata desirable properties as simplicity of -design, overall reduction in the effects of noise in 
the system, minimization of the order of the state variable estimator, and cost of implementation. 
For example, one possible method by which to invert the nonlinearities in any of the 
measurements 21001 to 21002^ is via a lookup table based upon offline calibration runs; another 
possible method, also based upon offline calibration, is via polynomial expansion. The latter is 
the method used in one embodiment of the present invention, as described in Detail 10 below. 
Noise reduction may be accomplished by filtering, for example by using finite impulse response 
(FIR) or infinite impulse response (IIR) digital filters, or else analog filters. The structure of an 
IIR noise reduction and data fusion filter, and its coefficient values may be determined by trial 
and error or by analysis. For example, a positional estimation filter could be designed via 
Kalman filtering techniques, in which a stochastic model of the input signal and state 
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measurement noise is combined with a model of the transconductance and transduction dynamics 
(such as equations (1 8)-(19) above) to resolve the order and coefficient values of the estimation 
filter. One skilled in the art will realize that various different filtering techniques can be used. 

Page 55, line 19, please amend the following paragraphs: 

As shown above, the nonlinearities in the electromechanical equations (18) and (19), 
which result from the position dependence of the L^ix), K(x) and Bl(x) , produce a nonlinear 
response in the transduction output jc as a functional of the voltage input V^^^,{t) . In-operation 
measurement of at least one position-indicator variable, together with suitable DSP computations 
as described above and in Detail 5 below, is used to calculate approximations to 
x{t) , x(t) , L^(x(t)) , K(x(t)) and Bl(x(t)) at any given moment during transducer operation. 
These numbers, together with the audio program input V^^^^^{t) , are then used by the controller 
circuit to implement a nonlinear feedback law for the transducer voltage input, F^oz/CO > based on 
the physical model of the system, as described by the control models given in equations (20), 
(22) and 0^ (40) . The overall control model obtained by combining the three control laws given 
by equations (20),_(22) and (3^(40), namely that given by equations (36) and (3^(39) above, 
was implemented in one embodiment of the present invention; the measured power spectrum 
distribution for a standard two-tone test, both with the combined correction and with no 
correction at all, are presented for this embodiment in Detail 14 below. It is seen that the effect of 
this combined feedback law is to eliminate or greatly reduce the distortions of the 3" Audax 
speaker transducer for which the data of Detail 14 were taken. Both intermodulation and 
harmonics peaks were significantly reduced. 

In the course of the derivation of the control laws in this section, it was noted that the 
physical audio transducer parameters L^{x),-K{x) and Bl{x) , as well as the position state 
variable x , are not perfectly known, and that for that reason, full correction as it appears in the 
equation of this section, will not in fact occur. The equations were derived assuming perfect 
knowledge by the controller; this was done to make the derivation of the control laws more 
transparent. In practice, however, these physical parameters and state variables are close 
estimates of their actual values. The attendant errors in modeling and measurement - both 
systematic and noise errors - introduce a small amount of unmod e led dynamics uncertaintv in the 
system. 

It is a well-known result in control theory that under certain conditions, unmodeled 
dynamics can lead to instabilities in a dynamical system under feedback. Care has been taken in 
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the implementation of the feedback laws of this section to reduce the sensitivity of the 
electromechanical system to these unmodeled dynamic s this uncertainty , thus preventing the 
possibility of dynamic instability in the electromechanical system, provided the coil/diaphragm 
excursion is not too high. 

Anyone skilled in the art will realize that other processes and process-components can be 
included in the transducer physical model, in addition to the transconductance and transduction 
which that are respectively encoded in the electric and mechanical equations (18)-(19). Examples 
of such additional processes are frequency partitioning and sound conditioning. These can be 
included in both the physical and control models, in accord with the modular approach to control 
modeling and implementation described in Detail 1 above. Similarly, the control models herein 
described can also be improved by accounting for ethe Fsmaller effects and terms within the 
electromechanical physical model, such as the terms that are not present in equations (18)-(19) 
but are present in equations (6) through (16). 

DETAILED DESCRIPTION 3: JUSTIFICATION OF APPROXIMATIONS 

A simplified physical model of a general speaker transducer, together with a modular 
collection of control models designed to implement linearization filters for sub-processes within 
the physical model, were presented in Detail 2 above. There are two ways in which these 
mathematical models are used in the context of the present invention: in actual physical 
implementation, and in simulation. 



Page 58, line 1, please amend the following paragraphs: 

As explained in Detail 1 above, there are many nonlinearities in the physical processes 
governing transducer operation, such as nonlinear elastic restoring force (i.e. nonlinear effective 
spring "constant"); nonlinear motor factor; nonlinear effective voice coil inductance; and motor 
BEMF, to name the most important ones. Computer simulations based upon the transducer-plus- 
controller model (and thus incorporating the leading nonlinear processes listed above) were used 
in the present work to study the effect of all of these nonlinearities, thereby elucidating the merits 
of implementing partial correction for a subset of the nonlinearities. For instance, it was found 
via simulation that transconductance nonlinearities (BEMF and inductive) are responsible for 
significant audio distortions at various important frequency ranges, which- led to the inclusion of 
corrections for these effects in -the control law (equations (20) and (22) above). In fact, 
dependent on program material, correcting for non-linear spring effects can have the 
consequence of increasing the excursions of the transducer coil/diaphragm assembly and thus 
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increase the non-linear effects of BEMF and 4 (x) . Nevertheless, it is still possible to achieve 

improved audio performance, especially at he low end of the audio spectrum, by correcting only 
for the nonlinearities in effective spring stiffness and in the motor factor. Thisy fact, as well, had 
been predicted by simulations of the model, and corroborated by experiment. 

We present several fevr^key simulation results relevant to the invention herein disclosed. 

FIG. 14 shows curves 4100 of simulated Power Spectral Density (PSD) which illustrate 
the effect of the transduction corrections alone (spring stiffhess and motor factor correction, 
equation (34)) both with and without BEMF and nonlinear inductance in the system. In FIG. 14 
the vertical axis is a measure of PSD in relative dB units. The curves of FIG. 14 were generated 
by simulating the performance of a particular transducer (that of the Labtec Spin 70 speaker) 
using a single 100 Hz tone; each curve clearly she¥m -shows that the highest power is in the 
fundamental 100 Hz tone, but that significant power is also present in the various harmonics of 
this tone. Overall, the curves of FIG. 14 shows that even at frequencies where BEMF is 
significant, introduction of corrections for spring and force constant greatly improve the system 
performance. Curve 4103 depicts the simulated PSD with no BEMF voltage term modelled, with 
linear (i.e. position independent) inductive EMF voltage term modelled, and with no correction 
incorporated in the modelling; the harmonics, and power present at non-harmonic frequencies, 
are an artifact of the finite time windowing used to perform the FFT (Fast Fourier Transform) in 
the simulation. Curve 4101 shows the PSD when the position-dependent (nonlinear) BEMF and 
position-dependent inductive EMF voltage terms both modelled, but still with no correction; the 
harmonics, as well as the general diffuse high-frequency content of the power spectrum, are seen 
to be enhanced by nonlinearity-caused distortions. Curve 4102, again depicting the PSD with 
nonlinear BEMF and nonlinear inductive EMF, but this time with transduction corrections, 
shows a marked decrease in harmonics and other, diffuse high-power spectral content. Finally, 
curve 4104 depicts the PSD with no BEMF and with linear inductive EMF, as in curve 4103, but 
with the difference that the transduction correction is applied. 

It is inevitable that there will be some delay between measuring and reading the sensor 
output, and sending out the command to compensate for the position-dependent nonlinear spring 
stiffness and motor factor (and for any other nonlinearities for which terms are included in the 
controller). Using model-based simulation, it was possible to determine that the existence of this 
delay, while somewhat degrading- the performance of the control algorithm, did not cause a 
significant problem, nor did it render the algorithm ineffective. 

The curves of FIG. 15 illustrate the reduction in distortion as a function of the delay in 
the correction loop 4200. As in FIG. 14, the vertical axis is a measure of relative PSD 
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magnitudes in dB. The curves of FIG. 15 depict the simulated PSD of the transducer-cone 
velocity, again for a 100 i/z audio input tone. In obtaining these simulation results, it was 
important to keep the amount of the nonlinearities the same for all the cases that were 
considered. This was achieved by suitably scaling the driving force as the time delay was varied. 
It is clear, from the curves of FIG. 15, that longer delays in the correction loop will increase 
distortion. However, for a 100 i/z tone, even at 200 //sec delay, the distortion is seen to be less 
than that of the uncorrected system. Curve 4201 depicts the PSD with no correction; curve 4202 
depicts the PSD with/ transduction) - correction but for the ideal case of no delay; while curves 
4203 and 4204 show the PSD curves with correction modelled and with simulated delays in the 
amounts of 100 //sec and 200 //sec , respectively. 

While a complete nonlinear spring cancellation will reduce the distortion in a_speaker's 
acoustic output, it will also remove the restoring force that was provided by the mechanical 
spring in the uncorrected speaker transducer, as discussed in Detail 2 above. In order to keep the 
speaker cone centered near its equilibrium position and place the mechanical resonance of the 
speaker at the desirable frequency, linear stiffness can be added electronically, as seen in Detail 2 
above. FIG. 16 displays a plot 4300 depicting the position of cone (i.e. the axial position of the 
coil/diaphragm assembly) in the presence of a single-tone excitation. Without the added 
electronic contribution to the effective spring stiflfiiess, the cone may drift from its equilibrium 
position, and may reach its limit of excursion; this is illustrated in the simulation shown in curve 
4302. Curve 4301 shows the corresponding simulated time-dependent cone excursion when an 
eIectronically-_added linear spring constant (suspension stiffness) is incorporated in the model. 

Page 61, line 18, please amend the following paragraph: 

FIG. 19 is a simulated version of spectral plot results 4600 of the two-tone 
intermodulation and harmonic distortion test for which actual, physical implementation results 
are reported in Detail 14 below. The two input tones are at 60 Hz and 3 kHz , and the portion of 
the simulated power spectrum distributions (PSDs) shown in the curves of FIG. 1 9 are in the 
vicinity of the- 3 ^/fe . The curves (4601 through 4603) clearly show the forest of intermodulation 
peaks, spaced uniformly 60 Hz apart and with decreasing power level away from the 3 kHz 
main peak. As is the case for the real spectrum in this frequency region (FIG. 65), the simulation 
shows the intermodulation peaks to be significantly suppressed when all four linearizing-filter 
corrections are applied (i.e. with the combined correction law given by equations (36)-(3^(39)). 
But unlike in the physical implementation, it is possible to select arbitrary time delays in the 
simulation. Two different delay values were chosen for this simulation: 10//sec and 50//sec . 
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And delays were only applied for the corrected runs. Curve 4601 shows the simulated 
uncorrected PSD; curve 4602 shows the dramatic intermodulation reduction when the 
corrections are applied, with 10 //sec simulated delay. Finally, curve 4603 shows the simulated 
PSD with corrections and with the longer simulated delay of 50 //sec . 

Page 62, line 10, please amend the following paragraphs: 

It will be clear to those skilled in the art that simulation of any particular implementation 
of the linearization and control methods described in this disclosure provides valuable 
information for practically implementing such systems for any particular application; and, 
furthermore, that the simulations developed here can be greatly expanded to cover many such 
systems and applications. 

DETAILED DESCRIPTION 4: STATE MEASUREMENT THEORY 

The present invention is described in the context of controlling an audio reproduction 
system, in part, by a model requiring real time measurement of at least one position-dependent 
state variable of the speaker transducer. In particular, one such state variable is the axial position 
X of the coil/diaphragm assembly. Real-time values of the state variable x are needed during 
transducer operation in order to effect the linearization of the transconductance and transduction 
processes, as set out in Detail 2. According to the present invention, it is unnecessary to have a 
direct measurement of jc ; it suffices to measure, instead, a position-indicator state variable, i.e. a 
variable which varies monotonically (but, in general, nonlinearly) with x within the range of 
possible diaphragm excursions. Once this position-indicator nonlinear state variable / (x) is 
calibrated against jc , real time measurements of the state variable / (x) can be used by the 
controller to effect linearization. 

Page 63, line 16, please amend the following paragraphs: 

The IR method is fully described in Details 8 and 13. The method is fiilly described in 
Details 6 and 11. The C method is fully described in Details 7 and 12. The position information 
derived by and C methods is generated using internal electronic parameters of the 

transducer. In contrast, the IR method is based on an external measurement of position. In all 
cases, to be useful as stand-alone position indicators;, the respective variables must be 
monotonic, but not necessarily linear, with position. It will be appreciated that there are other 
possible position indicators according to the present invention, which are measurable from 
internal electronic circuit parameters of the transducer that are not constant during transducer 
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operation, but instead vary monotonically with jc . One of ordinary skill in the art will readily 
recognize that there are many measurements that can be made on an audio transducer, but that 
K{x), Bl(x) , and L^(x) are commonly presented as the parameters most responsible for the 
noniinearities in the operation of such a transducer. The relationship of these parameters to these 
nonlinearities was explained in detail in previous sections, as was the fact that L^(x) also varies 

somewhat with frequency and depends on temperatures in the coil and within the magnet 
assembly. 

As an example of the use of position-indicator measurements in the controller in the 
context of the present invention, we consider one of the sub-process linearization laws presented 
in Detail 2 above; namely, the transduction-process control equation (34), where the transduction 
parameters S and B are non-constant functions of x . Any nonlinear position-indicator state 
variable / (x) can be substituted for x , as long as the positional related information is 
monotonic with x and is well behaved over the range of interest, i.e. the range of coil/diaphragm 
excursions in actual audio operation over which the correction is required. In other words, a 
nonlinear expansion in x can be replaced by a nonlinear expansion in any measurable variable 
that has a monotonic relationship with x over a suitable range of values. Thus, the variables S 
and B can be redefined as functions of x^^ , , or C^^^^„^,.^ , depending on the positional- 
detection method selected. The control law (34) then assumes the following different forms: 



i{t)^S(x,,)iwB{xJ 



mi{t) = S,ix,) + wB,{x, 




m 



i{t)^S{LjiwB{L^) 



i{i)=S{zJiwB{Z^ 



(41a) 



m 



(42a} 
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Page 65, line 11, please amend the following paragraphs: 

In accordance with the present invention, the transduction control law, equation (34), has 
been used to illustrate the use of nonlinear position indicators for linearization corrections. 
However, the same indicators can be used for some of the other corrections that can be added in 
a modular fashion to any particular implementation. These combinations of the modular control 
laws, described in the context of the present invention, are given by the control equations (20), 
(22), and (36)-(^{39} in Detail 2 above. In the case of the BEMF correction (equation (20)), the 
motor factor Bl(x) can be stored in the controller as a function of the nonlinear state variable 
f(x), while the instantaneous velocity x can be obtained not by measuring a motional state 
variable, but rather via numerical differentiation of the position, which in turn is obtained from 
/ (x) via the stored inverse functional relation /"^ . All controller-stored flmctions, whether 
having the form of polynomials, look-up tables or splines, or some combination of the these, will 
be computed, based upon calibration or characterization of the transducer, 'offline'; i.e. before 
actual transducer operation. 

Page 66, line 1, please amend the following paragraphs: 

Information from other external measurement apparata not utilized in the context of this 
invention, such as accelerometers, microphones, voltages from additional coils and/or additional 
transducers, can also be used to provide additional state variables, and thus can be used to add 
precision to, or reduce the noise, for positional or motional estimates. 

DETAILED DESCRIPTION 5: S AND B MEASUREMENT THEORY 

The present invention is described in the context of extracting the positional state of the 
speaker transducer's coil/diaphragm assembly, in operation, using measured state variables, from 
either internal circuit parameters, or signal(s) from external position-sensitive device(s), that are 
variables with that position. Measurement of all the parameters required to estimate S and B 
(the transduction-process variables introduced in Detail 2 above) with commercially available 
test equipment is both time consuming and fruitless. For a viable control scheme, the parameters 
must be regularly updated as they are sensitive to both time and temperature changes. 

Accordingly^ a method to measure S and 5 in a timely manner is described. The method 
used in this embodiment of the invention, and described in this section, to make the current value 
of B available to the controller DSP during operation, is also utilized for the electrodynamical 
transducer parameters Bl and Z^, as described in Detail 10 below. The values of Bl and are 
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needed by the controller in order to implement the transconductance corrections, namely the 
BEMF and inductance corrections respectively, as explained in Detail 2 above. 

Page 67, line 5, please amend the following paragraphs: 

Values for S can be measured directly from the control loop 6100. Considering the 
linearization correction equation (34) (or its subtracted version, equation (^£40}) for the 
transduction process alone with no audio information w, and hence without the B term, the 
spring force term S can be output independently simply by outputting a DC value - because for 
a DC signal, the only force in the correction equation is the static (spring-force to motor-factor 
ratio) term S(x) , and the numerical value of S can thus be measured. And since the 
corresponding numerical DC value of the arbitrary measure of position f{x) is also measured 
and fed back to the controller 6101, the approximate functional dependence of S upon / can be 
extracted via a suitable polynomial fit, and then used by the digital controller 6101 to look up the 
value of S which goes into real-time linearization correction of an actual, AC audio signal. 

FIG 2 1 is a flow diagram of a process for determining 5 as a function of position of a 
transducer. FIG 22 shows the voltage waveform 6206, the current from which is utilized to move 
the cone of transducer 61 03 and thus to determine and plot 5 as a function of x . Waveform 
6206 is output in step 6201 vr^rieh -that moves the diaphragm through positive and negative 
values of position jc , relative to the no-drive equilibrium value :c = 0 , over the range of the 
transducer's excursion. If, as is the case in the current embodiment-, a voltage-controlled 
amplifier -is used, a voltage ramp 6206 is output from controller 6101 as shown. After a new 
discrete voltage level on the ramp is output in process step 6201, a short wait for settling is made 
(process step 6202). The corresponding position-indicator state variable. /(jc) is then measured 
in process step 6203. The next discrete voltage level is then output in step 6201, unless a 'last 
step' decision is made in process step 6204; in which case the process ends with step 6205. Since 
a particular staircase signal is provided which is converted into the drive voltage V , and / (x) is 

measured simultaneously, this in effect constitutes the outputting of S(f'\f(x))) , i.e. the 
functional dependence S o of S upon /, where the circle symbol indicates function 
composition. The numerical value of the control parameter S used in the control loop 6100 is 
the transducer-coil current in voltage units - which is taken to be K . This procedure is 
approximately correct (in the case of a voltage controlled amplifier assumed here) to the extent 
that the non-Ohmic EMF terms in the coil circuit, including the effective coil inductance and 
BEMF voltage terms, are neglected. This is a justifiable approximation for sufficiently slow 
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ramping, i.e. long ramp-times and settling times. The ramp is made slow relative to audio signal 
timescales, because it is undesirable to put out audio information in the ramp. Therefore, the 
current into the coil is proportional to voltage by Ohm's law, to a good approximation. 

However, care must be taken that the ramp not be too slow, for otherwise significant 
heating of the coil could take place, and the coil current through the coil would then drop due to 
increased coil resistance. Care must also be taken to minimize the thermal and viscoelastic 
hysteresis effects reflected in the staircase-ramping measurements. Additionally, what 
unavoidable hysteretic effects do remain^ should be compensated for via some averaging 
procedure. In preparing the curve of 5 as a function of x for an Audax 3" transducer, waveform 
6206 shown in FIG. 22 included thirty-two steps of equal duration per each sweep from highest 
to lowest or lowest to highest voltage value. During the first and last of the steps the output 
voltage was zero. In each of the other steps, the voltage increment or decrement was 1 / 1 6 //? of 
the zero-to-peak amplitude of the waveform 6206, which was 0.25 volt. This value was before 
amplification. The amplitude of the ramp-sweep voltage signal fed to the voice coil of transducer 
6103 was about 20 times higher. This amplitude is determined, for each speaker transducer, by 
the need to cover the full excursion of the coil/diaphragm motion that is encountered in normal 
operation. 

Page 69, line 29, please amend the following paragraphs: 

As for the B term in the control equation G% )(40\ measurement of the functional 
dependence of B(x) upon / (x) , denoted as 5 o (/ (x)) , can be made by outputting a low 
amplitude tone, at a frequency sufficiently removed from the mechanical resonance frequency of 
the transducer to simplify the transducer's linear-response transfer function. The sound pressure 
output, or SPL, is measured at some fixed distance in front of the speaker, for example by means 
of a microphone, or alternately via other transducers within the speaker enclosure, or transducers 
in other speaker enclosures within a suitable proximity to the transducer being characterized. The 
off-resonance choice of tone frequency provides a relatively simple relation between the 
measured SPL and the motor factor Bl , which in turn is inversely related to B . The deduced 
values of B can then be tabulated against corresponding measurements of / (x) , for a stairway- 
ramped voltage signal 6206, in a manner similar to that used in the S measurements described 
above. At each DC voltage level, the low-l ew-amplitude tone is applied after that DC level has 
been held a sufficient time to allow electromechanical relaxation of the transducer to a steady 
state current and mechanical equilibrium. The frequency of the tone is fixed for eac hstairway 
stairway- ramped voltage sweep, but can be varied from sweep to sweep. However, the foregoing 
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approach is complicated by two factors. Firstly, the speaker's acoustic transfer function 
(diaphragm motion to SPL) is not a priori known for realistic speaker enclosures; and secondly, 
the suspension stiffness still affects the conversion of SPL values to B values, through the - 

dependent elastic resonance frequency, for tone frequencies low enough so that coil-inductance 
effects do not spoil the simple Ohmic conversion of voltage to coil current. This latter fact means 
that the S and B measurements are effectively entangled, as the extraction of B values requires 
knowledge of S values; and the converse also holds, as explained below. 

Because of these complications a hybrid approach is utilized, as follows. First, a Klippel 
GMBH laser-based metrology system is used to find an eighth-order polynomial fit to the 
function Bl (jc) , and the ratio function 

BM = ^ (44) 

where x 0 is the equilibrium position, is computed and replaced with a suitable lower-order 
polynomial fit. Note that this initial stage need only be performed once per given speaker, since 
drifts in the motor-factor function Bl {%) are almost entirely multiplicative, stemming from 
temperature dependence of the airgap magnetic field, and thus hardly affect the ratio B{x) . Next, 
a stairway-ramped voltage sweep of the type described above is performed, in which the 
position-indicator nonlinear state variable / {x) and the actual position jc are simultaneously 
measured. The latter is measured via a position s e nsor us e d with a Klipp e l GMBH m e trology 
system. This Klippel-tvpe laser, which returns a voltage known to vary linearly with actual 
position to a high accuracy. And finally, the Klippel-derived polynomial fit to B{x) is combined 
with the interpolated function / (jc) to yield an approximate polynomial interpolation for the 
composite functional relation B o f~\f{x)) : 

B o r' {f{x)) =b,^b, fix) + b, /{xf + b, fixf + . . . (45) 

Once interpolative approximations (polynomial or other) to both the functional relations 
S o/~* and B o (i.e. both S(x) and B{x) as functions of f{x)) are determined, these 
interpolations are stored and integrated into the controller DSP and used, in transducer operation, 
to dynamically compute and output a corrected coil voltage V^^.j from the original audio input 
signal w, via the control equation (3^(40), as explained in Detail 10 below. 

FIG. 23 is a general block diagram of a system 6300 depicting an audio transducer 6304 
with the digital controller 6301. Digital controller 6301 received two inputs: the audio voltage 
signal w 6302 (also referred to as ; see Detail 2), and the most recent measurement of the 
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position-indicator nonlinear state variable^ f{x) 6303. This nonlinear state variable is measured 
in the transducer 6304. Digital controller 6301 combines the audio input with the measured value 
of / {x) to compute the corrected V^^.^ in accordance to the control law. The control law may be 
that given by equation 0 ^(40) in the event that only the transduction-process corrections are 
selected, or by other equations in Detail 2 in case the user decides to activate other combinations 
of control laws. The voltage F^^,, is output in analog form 6305 by digital controller 6301, and 

provided to the amplifier 6306. The output voltage fi-om amplifier 6306 is provided to transducer 
6304. 

As discussed in Detail 2, the use of the entire spring force in the correction, thus in effect 
electronically subtracting away the entire elastic restoring force, would lead to dynamical 
instability. It is therefore necessary to add back a linear spring restoring force calculated as an 
adjustable fraction of the measured spring factor at equilibrium, 5(0) . This is done by 

subtracting a term linear in the estimated position (/ {x)) from the ratio of the S o /"^ (/ (jc)) 
polynomial to the B o f~\f{x)) polynomial, since this ratio is the constant times an 
interpolating function for the suspension stiffhess x K{x) . The net result of this subtraction is 
that the numerical values of S , and the functional relation S o , are replaced by new 
quantities, denoted here as 5" and 5" o respectively, in the control equation (3^£39}. If the 
transconductance corrections are turned off, equations (36) and 07^ (39) reduce to the 
transduction-corrections equation {^{40}, which is just equation (34), but with S replaced with 

the following subtracted value: 

S'=S-kr\f{x))B 

— Wbere -wherca -A: = q K{Q) I is a constant multiplier^ related to the adjustable parameter 
q of equations (37) and (3 8 ). (39) and (40). The multiplier q can be optimized by user 
preference. In Equation (46), the three quantities S , B and S' are all expressed as interpolated 
polynomials in the measured position- indicator nonlinear state variable / (jc) , as described 
above. 

Page 73, line 9, please amend the following paragraphs: 

According to the present invention, there are provided parameterized linearization-filter 
functions characterizing the given transducer, which are measured and estimated using in- 
operation measurements of at least one nonlinear position-indicator state variable, augmented by 
preliminary (characterization) calibration runs in which this nonlinear state variable is measured 
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simultaneously with a more linear position-indicating variable (such as the Klippel-GMBH laser 
metrology system). The nonlinear position-indicator variable measured in operation can be a 
voltage output from an optical device, as is the case in one embodiment of the present invention 
and as is described in Details 8 and 13 below; or it could be an output from the internal 
electronic parameter measurements, as described in Details 6,7,1 1 and 12. These measurements 
could be augmented by an external measurement of sound pressure level during characterization 
runs, as described above. 

Accordingly^ an invention where t he S and B parameters, which are needed by the 
controller to implement the transduction-process portion of the linearizing control law, can be 
matched to the program material by adjusting the parameter q governing the electronic spring 
force compensation, as described in equations (37V (3 8 ) (39). (40) and (46). 

DETAILED DESCRIPTION 6: MEASUREMENT THEORY. 

An important aspect of the present invention is described in the context of a digital 
control system which linearizes audio reproduction using a position-indicator state variable, f(x), 
which is monotonic in position. The inductance of a transducer voice coil provides such a 
position -indicator state variable. This method appli e s to many other class e s of non lin e ar 
actuators and motors. 

Page 74, line 11, please amend the following paragraphs: 

It can be seen from curve 403 in FIG. 4 that the values of (in this case at frequencies 
below I kHz ) are monotonic with x ; that is to say, no two distinct x values within the range 
- 2 mm to 2 mm correspond to the same value of . We can thus map Bl (curve 401 ) and K 
(curve 402) onto 4 » ^'^d a measurement of will uniquely predict both Bl(L^) and K{LJ. 
These functional relationships are depicted in FIG. 24, in which curve 5101 is a plot of K in 
Newtons/mm and curve 5102 is a plot of Bl in Newtons/amp, both of which are plotted against 
4 for the same data from FIG. 4. This new mapping provides a basis of a correction scheme. 
Because the inductance of the voice coil is a function of its position, by measuring the 
inductance the position of the voice coil is determined. Thus provides an inductive position 
detector. 

Page 77, line 9, please amend the following paragraphs: 
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A method for measuring the coil inductance is illustrated by the block diagram in FIG. 
31. A supersonic probe tone ("carrier signal") is applied via input line 7401 to the voice coil of 
transducer 7402. In this approach, a reference R L circuit 7403 is placed in series with the voice 
coil. The supersonic signal is then injected into the voice coil of the transducer 7402 in addition 
to the audio signal, and the voltage across the voice coil of the transducer 7402 and the reference 
R L circuit 7403 is measured. Reference R L circuit 7403 may be implemented using a resistor 
and a coil in series. Alternatively, a coil or a resistor may be used to implement circuit 7403. The 
measured voltage signals are sent via summer summing junction 7404 and s umm e r summing 
junction 7405 through filter 7406 and filter 7407, respectively, and the ratio of the output of the 
filters is then determined in either the analog or digital domain. The filter 7406 and filter 7407 
are band pass filters implemented about the fi-equency of the carrier signal. Envelope detection 
via envelope detector 7408 and envelope detector 7409 is used to extract the signal due to 
changes in . The ratio of the vohages coming out of the envelope detector 7408 and detector 
7409 can be described in the Laplace domain as: 

y~.-Y^ (47) 

— Whef ewhere /?/ is the resistive component of coil impedance at the probe tone frequency, 
including both the Ohmic coil resistance R^ and the lossy effective coil impedance component 
due to eddy currents. R^j- and Z^^ are the respective series resistance and inductance of the 
reference R L circuit 7403; and s is the Laplace variable. Because the ratio of the two voltages 
is taken, the signals that are close in frequency to that of the carrier, and thus cannot be rejected 
by the band-pass filter 7406 and filter 7407, will not introduce significant error in 

L R 

determination. As long as L^^^ and R^j^ are chosen so that — — and — — are the same for 

^ref ^ref 

R L 

frequencies near the probe tone, V^^^^^ remains a constant equal to — — = — — , regardless of the 

Kef ^ref 

presence of -other signals in the system that are close to the frequency of the carrier signal. Since 
4 varies with coil position x, V^^^.^ will change accordingly. FIG. 27 shows the Bode plot of 
the transfer fiinction F^^^.^ given in equation (47), while FIG. 28 shows the corresponding phase 
Bode plots. 

The ordinate in FIG. 27 is the magnitude of K^^,.^ , in dB units, while the ordinate in FIG. 
28 is the phase of F^,^^ , in degrees; in both plots, the abscissa represents angular fi-equency in 
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units of radians per second. In both FIG. 27 and FIG. 28, the family of Bode plots is for 
progressively larger values of 4 , with the highest 4 value resulting in the curve 7201 and 
curve 7204, while the lowest value results in the curves 7202 and 7205. It is seen that as 4 
increases, so does the magnitude of F,^,,^ . The sensitivity of F^^,.^ to changes in 4 is clearly a 
function of the probe tone frequency. The higher this frequency, the more sensitive V^^^.^ will be 
to 4 variations. 

To reduce the effect of the common mode in-band noise, which is present in the voltage 
across the voice coil (i.e. (4 ^ + * O and the-in the voltage across the reference R L circuit 

(i.e. (4e/ s -f R^^.) • / ), upon the voltage ratio, the phase shift of ^ must be small. 

Thus, the choice of probe tone frequency may have an impact on the effectiveness of noise 
cancellation within the above-described approach. Furthermore, to iftsuf eensure the noise 
cancellation advantage of this algorithm, the band pass filters, mentioned above, must be 
matched as closely as possible. 

Page 79, line 22, please amend the following paragraph: 

Because FIG. 27, FIG. 28, FIG. 29 and FIG. 30 illustrate that a thermal change in R^ is 
likely to be misinterpreted as a change in 4 , a modification in the algorithm is needed to 
separate this thermal effect from actual changes in 4 that are caused by changes in the voice 
coil position. From FIG. 29 and FIG. 30, it is clear that the effect of variations in R^ upon the 
ratio F^^^,^ is minimized at the higher probe tone frequencies. This characteristic of V^^^.^ can be 
utilized to accurately determine in the presence of thermal changes to R^ . For instance, for 
the Labtec Spin 70 speaker transducer for which the curves in FIG. 27, FIG. 28, FIG. 29 and 
FIG. 30 were generated, the use of a carrier signal at \5QkHz will significantly reduce the 
thermal effects upon measurement. 

DETAILED DESCRIPTION 7: C THEORY — PARASITIC CAPACITANCE 

AND CANT DYNAMICS. 

Page 80, line 8, please amend the following paragraphs: 

An important aspect of the present invention is described in the context of a digital 
control system whieh -that linearizes audio reproduction using a position-indicator state variable, 
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/(jc), which is monotonic in position. The parasitic capacitance C^„^^„^,^ between the voice coil 

and the body of a transducer can be used to give such a position -indicator state variable. This 
method applies to many other classes of non-linear actuators and motors. 

The parasitic capacitance C^^^^,,^^ between the voice coil of a transducer and the body of 

the transducer is largely determined by the relative positions of the voice coil and the magnetic 
pole pieces and central core. The variation of this capacitance with position is relatively 
straightforward and robust (reproducible). As illustrated, for example in FIG. 3, typically the 
voice coil 303 fits about a central core -310 which that is part of the iron assembly 305. The 
variation in the parasitic capacitance depends largely on the overlap of the voice coil 303 with 
the central core 310, and, to some extent, with the outer pole piece 3 1 1 as well. 

Page 81, line 4, please amend the following paragraph: 

Importantly for the purpose of the present C method, C^„^,,^^ is an easily measurable 

internal circuit parameter of the transducer which is, at the same time, a state variable which 
depends monotonically upon axial coil position x . As the coil moves deeper into the magnetic 
airgap, the capacitativ e capacitive contact areas between the metallic surfaces of coil and poles on 
the one hand, and between former and poles on the other hand, increases; and thus so does the 
value of the parasitic capacitance. 

Page 83, line 2, please amend the following paragraphs: 

It is possible to understand the results from FIG. 32 and FIG. 33 using simple semi- 
quantitiv e quantitative models. Although some fairly involved modeling is required to obtain an 
accurate prediction of C^^^^^,^,^ {x) for a given transducer, it is quite easy to estimate its order of 

magnitude. Thus, referring to FIG. 35, assume a coil of height h and radius r . In the modeling 
described below, it is assumed that the coil former is to be non-conducting; thus only the coil- 
slug contribution to the capacitance is considered. Furthermore, we ignore the capacitance 
between coil and outer pole (as the capacitive overlap area for that pair of conductors is assumed 
smaller than between coil and core). For simplicity the wire indentations and insulation are 
likewise ignored. The maximal value of C^^„,,^(jc) occurs when x is smallest, that is to say, 

when the coil is farthest into the magnetic airgap, x = x^-^ . Assuming that at this coil position 
the capacitative c apacitive contact area between coil and slug equals the total area of the coil's 
cylinder, the following estimate -results: 

^ parasitic 
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— Whef ewhere is the permittivity of air and gin^er/or an estimate of the average distance 
between the steel of the central pole, and the copper surface of a typical wire belonging to the 
coil's innermost winding layer. For instance, in the case of the LabTec Labtec speaker transducer 
discussed above, the geometrical parameters are estimated to be r = 1.5 mm , h - 5mm , and 

Sinterior ^ 0.2mm ^ SubstitutioH of these three values into equation (48) yields: 

C^rasi,c(Xn,n)^^OpF (49) 

The value measured electronically was found to be about ISpF for this transducer. The 
discrepancy is reasonable given the parameter estimates. 

For transducers of smaller speakers, such as those utilized in cell phone receivers^ smaller 
capacitance values, for example several picoFarads, were measured. This decreased magnitude 
can readily be understood from the way in which the right-hand side of equation (48) scales 
down with the linear dimensions o f the speaker's transducer. 

Page 84, line 18, please amend the following paragraph: 

When an aluminum former is used as a heat sink for the voice coil, which is often the 
case in transducers of woofer speakers due to the high power levels dissipated in their coils, 
unwanted circumferential eddy currents are induced in the former. These eddy currents result 
from two effects: one is the EMF induced in the former due to the its axial motion through the 
radial magnetic field in the airgap; and the other is the EMF induced by the time dependence -of 
the coil-current's contribution to the axial magnetic field through the former's interior. In order 
to suppress these eddy currents, it is standard practice to interrupt them by introducing a slot 
along the axial length of the former's surface. This practice does not, however, completely 
eliminate the former eddy currents, but instead has the effect of distributing them nonuniformly 
around the former's circumference. These nonuniform currents, in conjunction with the static 
radial magnetic field in the airgap, cause magnetic Lorentz forces on the coil/diaphragm 
assembly wh4eh -that lack azimuthal symmetry. These non-uniform forces lead to a non- 
vanishing torque, and therefore to canting. This former-caused canting effect is discussed in J. 
Vanderkooy, J. Audio Eng. Soc, Vol. 37, March 1989, pp.1 19-128. 

Page 85, line 14, please amend the following paragraphs: 

For the Labtec Spin70 transducer, canting due to fractional turns, in addition to 
exacerbating audio distortions, makes correction using the C method less desirable in some 
ranges of cone movement, by causing the function C ^^^^ {x) to become non-monotonic when 
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in operation. As the voice coil moves towards the back of the speaker through the airgap to, or 
beyond, its mechanical equilibrium point, the fractional wire-turns approach the region of high- 
magnetic-field in the airgap sufficiently to cause significant torque and canting; the cant, in turn, 

causes some parts of the coil wire's conducting surface to recede further from one or the other of 
the magnetic pole structures, increasing the value of the effective capacitive gap gi„tg„^^ in 
equation.(43) and thereby decreasing the values of C . 



Page 86, line 14, please amend the following paragraphs: 

It is assumed that the fractional part of the number of coil-wire windings is 1/2, and the 
above notation for coil dimensions is retained. A further simplification is made, in that the radial 
magnetic field at the position of the half-winding is replaced with the same field component 
averaged over all the coil's windings. The canting torque on the coil/diaphragm due to the 
magnetic Lorentz force, is then approximately: 

T„.^eac^J^Bl(x)iit) (50) 

— Wher ewhere r denotes torque; /(/) is the coil current, time independent in the DC case; 
Bl(x) is the transducer motor factor, and the total number of windings in the voice coil. 

This magnetic torque is opposed by an elastic torque, caused by the elastic restoring 
forces acting to counter the canting. We denote by jh^ PehsuM) ^(^) the relevant torsional 
spring constant - i.e. the elastic torque, per radian of tilt, exerted by the speaker's spider and 
surround upon the coil, diaphragm and cone; here h is the coil's height (defined above equation 
(48)), K{x) is the coil/diaphragm suspension stiffhess recognized in prior art, while (x) is 
a dimensionless elastic ratio modulus characteristic of the coil/diaphragm assembly. The p^i^^ti^ 
ratio modulus is expected to be significantly larger than unity, as speaker diaphragms are 
designed to resist canting while allowing axial motion. 

With the above definitions, the elastic restoring torque is simply: 

W« ^P,astM)K{x)e{t) (51) 

— where 6{t) represents the canting (or tilt) angle, in radian units, as a function of time. 

When the coil is driven with a DC or quasi-DC current, mechanical equilibrium is 
attained when the magnetic and elastic torques balance: this occurs at a tilt angle of 
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Ignoring the coil-wire insulation, this tilt results in an increase in the parasitic 
capacitance, roughly estimated at: 

1 J , 1^(01 



(53) 



— Wbef ewhere 1 9{t) \ is the absolute value of the tilt angle, and C^^,^^„,^ (jc) is the capacitance 

for the case of no canting. 

Since the driven-coil measurements for the LabT e c Labtec speaker transducer were 
quantified in terms of coil-circuit voltage rather than coil current, we set i{t) = V^^^i(t)/ in the 

above equations, where is the coil's Ohmic resistance (this relationship requires corrections 
in the AC case, as detailed elsewhere in this document-). Thus, for the DC case, equations (52)- 
(53) now yield the predicted fractional increase in- parasitic capacitance due to canting: 

^ ^ ^ y ^ 1 BI{X) ^^^^ 

C parasitic StT R^ N p^^^^^^ {x) K{x) 



Page 88, line 16, please amend the following paragraphs: 

This tilt angle would only result in a maximal lateral displacement of order 0.02 mm for 

parts of the coil - too small to cause the coil to be physically blocked by the pole structure, but 
enough to result in discernible audio distortions. However, the estimate for the fractional change 
in stray capacitance is quite dramatic, and in agreement with the measurements made for this 
speaker transducer. 



DETAILED DESCRIPTION 8: IR DIODE MEASUREMENT THEORY 



Page 90, line 7, please amend the following paragraph: 

The present invention is described in the context of controlling an audio transducer 
system in part byAsystem consisting of hardware and software. 



Page 90, line 12, please amend the following paragraphs: 

A DSP based controller 10101 consists of a PPS DSP processor and software system 
10102 and an interface system 10103 consisting of analog input/output and user interface 
software. Audio input is provided to DSP based controller 10101 through a signal-matching 
network 10104 w4rie hthat filters the audio input and provides the correct level of input to the 
interface system 10103. The audio input is acted on by the control routines in the DSP based 
controller 10101 and is output to a second signals-matching network 10105. The signal from the 
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signalrmatching network 10105 is provided to a power amplifier 10106. The output of power 
amplifier 10106 drives a speaker transducer 10107. A position sensor 10108, or sensors, is used 
to provide a position indication signal, indicating the position of the coil/diaphragm assembly of 
the speaker transducer 10107 to sensor signal conditioner 10109. Such position sensors could be, 
for example, the detector of Detail 6, or IR detector described in Details 8 and 13, or C 
detector described in Details 7 and 12. Sensor signal conditioning system 10109 is used to 
amplify and filter the positional signal and match it to the level required for the interface system 
10103. 

FIG. 39 is a block diagram of a particular embodiment of an audio reproduction system 
15100 that includes a DSP based controller 10101. A Personal Computer (PC) 15101, which 
could be an eMachines T1742, is used as a control and user input environment for the DSP based 
controller 10101. The DSP based controller 10101 is implemented using a M67 DSP board 
15102 and a A4D4 I/O board 15103 both manufactured by Innovative Integration Inc. (Simi 
Valley, CA). The M67 DSP board 1 5 1 02 is a mother-board for the A4D4 I/O board 1 5 1 03 . The 
M67 DSP board 15102 contains a 106MHz TMS320C6701 floating point DSP manufactured by 
Texas Instruments and has been modified to add an inverter (74LS14) between JP14 pin 34 to 
JP23 pin 29. The A4D4 I/O board 15103 consists of four 16 bit analog-to-digital converters 
(ADCs) and four 16 bit digital-to-analog converters (DACs) with interface circuitry to the M67 
DSP board 15102. A Lynx L22 card 15104 manufactured by Lynx Studio Technology, Inc 
(Newport Beach, CA) installed on the PC 15101 provides an audio signal 15105 whieh -that is 
input to the A4D4 I/O board 1 5 1 03 . The Lynx L22 card 1 5 1 04 receives input via Cool Edit Pro 
software 15106 (version 2) installed on PC 15101. The Cool Edit Pro software 15106 generates a 
'.wav' type digital sound file from a music source, which could be a CD player 15107 also 
installed on the PC 15101. After processing by the DSP based controller 10101, the corrected 
analog audio signal 15108 is output from the A4D4 I/O board 15103, and provided as an input to 
a 20:1 attenuator 15109. Output from the attenuator 15109 is provided as input to a Marchand 
PM224 amplifier 15110 with internal jumpers set to give a DC coupled amplifier. The Marchand 
PM224 amplifier 151 10 is manufactured by Marchand Electronics Inc (Webster NY). The 
Marchand PM224 amplifier 151 10 is used to drive a 3" transducer 15111 manufactured by 
Audax (Westlake Village, CA). The embodiment of audio reproduction system shown in FIG. 39 
uses the IR method of position sensing. An IR detector 151 12, the operation of which is 
described in Details 8 and 13, is used both to measure the position of the coil/diaphragm 
assembly of the 3" transducer 15111 and to match the signal to the input stage of the A4D4 I/O 
board 15103. The output 15113 ofthe IR detector 15112 is an input to the A4D4 I/O board 
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15103. 

DETAILED DESCRIPTION 10: SOFTWARE AND PROCESS FLOW 

Page 92, line 6, please amend the following paragraphs: 

FIG. 40 shows the process flow used to linearize the transconductance component of the 
signal conditioning process and the transduction process of a given audio transducer, based upon 
the control model given by equations (36)-(3^39) in Detail 2 above. FIG. 40 applies also for the 
case that only a subset of these correction s are applied. 

In the process illustrated by FIG. 40, the first step 1 1 1001 entails measuring large signal 
(LS) transducer parameters. This step yields coefficients of polynomial interpolations for the 
functions Bl{x) and L^{x) . The measurements are performed using a Klippel GMBH laser 
metrology system, with procedure as detailed in Klippel System Manual dated May 2, 2002. 

In a second step 1 1 1002, a software control program is invoked , for e xampl e the software 
control program in file 071 1 19.txt included in the computer program listing appendix provided 
by th e compact disks includ e d with this application . In a third step 1 11 003, the invoked software 
control program is run in 'Calibrate' mode in order to calibrate the fiinctional relation between 
coil/diaphragm position x and the position-indicator nonlinear state variable, /(jc), which in one 
embodiment of the present invention is the voltage output of the IR circuitry: x^^ = f(x). During 
this calibration, the software control program collects corresponding values of x as measured to 
an approximation by the Klippel laser, and / (x) , in relation to the corresponding values of 
voltage outputs^as described in Detail 5 so that the dependence of f{x) with x and the 
dependence of S with / (jc) can be determined. 

An example of the software control program used in step 1 1 1003 is provided by FIG. 41, 
FIG. 42, FIG. 43, and FIG. 44. The data obtained from steps 1 1 1001 and 1 1 1003 are used to 
find Best Fit coefficients for lowest order polynomials of *? , x , Bl and as fiinctions of jc.^ , 
as indicated by step 1 1 1004. Here 'Best Fit' is defined as that curve which is of the lowest order 
and which does not exceed specified rms and maximum errors, subject to substantial weighting 
in the mid section of the range of the / (x) variable. More details and specifics on 'Best Fit' are 
provided later in this section. The user then inserts the polynomial coefficients obtained from 
step 1 1 1004 into the Software Control Program - step 1 1 1005. Next, the user invokes the 
Software Control Program for Normal operation - step 1 1 1006 - and operates the program in 
Normal mode 111007. 
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FIG. 41 shows the structure of one embodiment of the Software Control Program that is 
used both for obtaining data during calibration 1 11 003, and for operating in normal mode 
1 11 007 in which linearized sound is produced. The initialization process II 1 101 places the 
system in a known state. The software control system can then be selected to operate in 
calibration mode 1 1 1103, which consists of an S and an x calibration process, or to operate in 
the normal mode 11 11 04. Typically, the first time around the user needs to select the calibration 
mode 1 1 1 103, as indicated in 1 1 1003. After completion of calibration mode 1 1 1 103, the system 
can be selected for normal operating mode 1 1 1 104, in which the software controls the sound 
reproduction process through an Interrupt Service Routine (ISR) 11 11 06. Note that the ISR 
fiinctionality 1 1 1 106 is also used in calibration mode. On an exit event 1 1 1 105 prompted by the 
user, the system stops the program 1 1 1 107. FIG. 45, FIG. 46, FIG. 47 and FIG. 48 cover the 
normal operation in detail, while FIG. 42, FIG. 43 and FIG. 44 cover the calibration mode in 
detail; all these figures are described later in this section. 

Page 97, line 5, please amend the following paragraphs: 

The operations within step 1 1 1305 are detailed in FIG. 51. In step 1 1 1401, for the given 
set of coefficients, for example, , c, , . . . Cg for a polynomial Y - values of Y^^^^ are calculated 

as follows: 

— Fof for each of several points p in the range given above. The above Y^^.^(p) values are then 
used in step 1 1 1403 to compute new coefficients, mi in module 1 1 1404 to compute errors. Here 
7^,^ values: 7^,^^^, Y^^^^^^ • • ^ong^^ are calculated for 33 points /?1, ;733 distributed 

uniformly over the above range. It will be readily recognized that that the number of points used 
can be changed within the fi*amework of this invention. 

Page 97, line 30, please amend the following paragraphs: 

Let Zj = Wj {Y^rig)j''> J = » be the weighted desired output vector, and b^^ , 
k = 0,...,M be the reduced order vector of coefficients \s4Heh -that needs to be determined. Then 
the weighted output vector for the points pj for the coefficient column vector b is given by the 
new column vector Ab . The total weighted squared error between the two weighted vectors is 
given by: 

E = (z-Abf (z-Ab) 
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Page 99, line 18, please amend the following paragraphs: 

The Matlab program dev e lop e d by Tymphany Corporation to implement modul e 1 1 1 103, 
is named 'reduce_order_of_XlsrSB_ 

Lcoeffs'. The code is included in the computer program listing appendix on the compact disks 
included with this application in files: 071 1 15A.txt; 071 1 1 15B.t>a; 071 1 11 5C.txt; 071 1 1 15D.txt; 
071 1 nSEixt; 071 111 5F.txt; and 071 1 115G.txt File 071 1 15A.txt li s t s the main program itself, 
which calls the filenaming utility function 'fNmsInOutXlsrSB Lcoeffs' (listed in file 
071 1 1 1 5B.txt) for the purpose of letting user specify input and output filenames. Next, the main 
program calls th e function 'r e duce ord e r of coeffs', (fil e 0711 1 1 5C.txt) for e ach of the 4 curves 
whos e order is to b e r e duc e d. This function in turn calls ' WtdLstSqPolyCoeffs' (file 
071 1 1 15D.txt), a subfunction that calculat e s th e coeffici e nts according to th e weighted l e ast 
squar e s approach d es crib e d abov e in ord e r to provid e th e best fit co e fficients. In turn, the 
subfiinction 'WtdLstSqPolyCooffG' of file 071 111 5D.txt call s 'cnstrct_A_^mat_a_2_cor (listed in 
fil e 071 1 1 15E.txt) to construct th e weight e d - y t matrix and th e w e ighted z column v e ctor n ee d e d 
for calculating th e co e ffici e nts. Th e function 'r e duc e ord e r of coeffs' of fil e 071 1 1 15C.txt also 
calls a plotting routin e , 'plt__data_sup_y3' (fil e 071 1 1 15F.txt), which plots curv e s if request e d by 
us e r. Finally, th e program calls th e function Svrls e tRdcdCo e ffToOp e nFil e ' (list e d in fil e 
071 1 1 15G.txt) four tim e s, in ord e r to writ e on e set of co e ffici e nts for e ach of th e four curves. 

Page 101, line 8, please amend the following paragraph: 

FIG. 47 shows a flow diagram of the ISR 1 1303. When the negative edge of the sampling 
clock occurs, the software control passes from the wait loop and command parser 1 1204 to step 
1 1401 . Step 1 1401 limits the value of the word to be sent to the DAC 1 1402, so that it does not 
exceed the input range of the DAC 1 1402. The DAC can be an onboard DAC as it is with the 
Innovative Integration A4D4, or a serialr-port-based off-board DAC. The analog signal that is 
created is the corrected audio signal V^^.j , and is fed to a power amplifier 10106. To create the 
corrected audio sample, the ISR module 1 1303 uses IR sensor data f(x) from module 1 1403 
and audio data from module 1 1405. A digital filter 1 1404 is used to minimize sensor noise in the 
measurement of f{x) . Module 1 1406 computes 5 , 5 , and corrections from the filtered 
value of / (jc) 1 1404, as described below. 

Page 102, line 3, please amend the following paragraphs: 
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The corrected audio signal V^^., , calculated by a combination of actions by modules 
1 1406, 1 1407 and 1 1408, is derived from input audio signal and the value of filtered f(x) using 
the following eight equations: 

Bl = Bl, + Bl, fix) + Bl, {f{x)f -f Bl, (fix))' (5357) 

S = S,+S, fix) + S, ifix)f + ... + 53 ifix)y - k ifix)) I Bl (#^58) 

k=L,+ L, fix) + L, ifix))' + L, ifix))' (5560) 
m = akt-r) + fiif-\fix, (t))) -f^' ifix^ it - r))) ) (5661) 
DEMF - iK^, Bl - IBl) i(/) 5£MF = --iKy^ I Bl) kit) iS^62) 
I r,,,,,(0/^/ I jygM/^ V,it)^S^V^^,,St)BhlBl + BEMF 

(5«63) 

Vco, it) = V, it) ^K,,L, iV, it) -V,it- T)) + K,, Bl m (^64) 
— where: V^^^j is the corrected voltage signal applied across the voice coil and including all four 
corrections (S, B , BEMF and ); is the corrected voltage without the inductive correction; 
^audio is the audio input voltage signal, suitably normalized; / and r denote the current time-step 
and the sampling time, respectively; and the constant k in the subtraction term in the polynomial 
expansion for 5* (last term onright-hand side of equation (#^{58}) is the electronic linear spring 
stiffness remaining after the linearizing filter (see Details 2 and 5 above). It is used in the 
calculation of S in order to maintain an appropriate level of restoring force in a transducer (see 
Detail 5 above); without this restoring term, the transducer would become unstable. 

Equation (§4){59} is a correction applied to linearize the IR position-indicator state 
variable x^^ = fix) if necessary. Equation {55)£60} is the correction for nonlinear inductance . 

Equation (&€) (6l) is a digital filter designed to estimate the velocity of the transducer 
needed for the BEMF correction. Equation (57){62} calculates the required BEMF correction. 
The BEMF correction comprises two components: the removal of the nonlinear BEMF and the 
replacement with a linear BEMF. The equations incorporate a multiplier for each term to allow 
for fine adjustment of the correction. Equation ( ^(63) and (§9){64} implements the above 
components of the audio correction. 

It will be appreciated that there are many different ways of discretizing the numerical 
differentiation operation of the control diagrams FIG. 1 1 and FIG. 12, and that the 
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implementation of these numerical differentiations used in one embodiment of the invention, and 
shown in equations ^ ^(61) and ( ^(64) , represent but one possible choice. 

Digital filters may be added to equation {S9)(64} for smoothing, equalizing and noise 
reduction. The polynomial coefficients as well as the powers of filtered / (jc) are stored in 
arrays, so that the needed sum of products can be easily computed. Moreover, the array for 
powers of filtered f(x) may be constructed recursively, again reducing the computational cost. 

Page 103, line 25, please amend the following paragraphs: 

FIG. 42, FIG. 43 and FIG. 44 show flow diagrams of S and x versus / (x) calibration 
111103. 

For calibration, the mainline loop is finite (while that in- normal mode is infinite) and 
results in a tabulated output, from which a polynomial curve is fitted and polynomial coefficients 
extracted for use in the Normal Mode 1 1 1 104. 

Page 104, line 7, please amend the following paragraphs: 

FIG. 43 shows the details of HW and ISR operations for S calibration 1 1504. It depicts 
Sampling Clock 1 1601 and ADC Convert 1 1602, which are similar to corresponding modules in 
FIG. 4-1-346; the same description applies, and is thus not repeated. Modules 1 1604 and 1 1605 
limit and convert the digital values to analog waveform. Module 1 1606 tests whether the data is 
to be collected. During calibration mode, the mainline S calibration loop 1 1505, detailed below 
in FIG. 44, sets and clears the flag 'Collect_data'. If this flag is set, the data collection proc ee ds 
asjs^done by the module 1 1607, and a sample count is tallied. Also, module 1 1608 reads the S 
value from the array, to be used in the variable 'dacvalue'. If the flag is not set, these two 
modules are bypassed. Module 1 1609 executes the return from ISR. 

FIG. 44 shows the details of mainline S calibration loop 1 1505. Module 1 1701 checks 
whether any value of S is left with which to operate the loop. If there is one, it executes the path 
comprising modules 1 1702 through 1 1707 to send out the S value via the ISR 1 1603, and to 
collect the corresponding value of /(jc) and x as follows. Module 1 1702 executes a wait of 
100 milliseconds to allow the transients in the transducer to attenuate. Module 1 1703 sets the 
'Collect_data' flag which signals the ISR 1 1603 to collect data. Module 1 1704 allows 1 
millisecond to collect samples, which at 48 kSPS collects 48 samples. These samples suffice to 
give a good reading of / (x) , the IR data^^ and x, the laser data. Module 1 1 706 performs 
averaging, and the module 1 1 707 stores S , f{x) and x for offline curve fitting. As long as 
there is an S value to be covered, the process continues. 
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Page 105, line 17, please amend the following paragraphs: 

FIG. 48 illustrates the details of the Wait Loop and Command Parser 1 1204, shown in 
FIG. 45, which is abbreviated below as WLCP. The system enters into the 11 801 step of WLCP 
from Enable ISR Setup and Enable Sampling Clock 1 1202; in step 1 1801 it is determined 
whether Normal Mode operation should stop. If 'Yes', system enters into step 1 1803, in which 
Interrupt is disabled and the HW is put into a known state; then system is passed out of WLCP 
and into User Mode Select 1 1 1 102. But if the answer to the 'Stop?' query (step 1 1801) is 'No', 
the DSP passes to step 'Command?' 1 1802, in which the WLCP checks to see whether User has 
entered a keyboard command since the last check (checks are spaced several microseconds apart 
during the Wait Loop). If no new keyboard command has been entered during the most recent 
such time interval, this is interpreted as a 'No' response to the 'Command?' query-, and the 
system is looped back to this 'Command?' query 1 1802. But if and when WLCP finds that a new 
keyboard command has been entered during the most recent time interval, each of the following 
optional keyboard responses are interpreted by WLCP as a 'Yes' and acted upon. User keyboard 
response 'c' causes the DSP to begin implementing corrections: 'Corrected Audio Mode' 1 1804; 
after this mode is entered, the system is passed back to the 'Stop?' query 1 1801 . User keyboard 
response 'b' causes the DSP to enter the mode 'Adjust Linear BEMF' 1 1805, from which it is 
again returned to 'Stop?' query 1 1801. The following are the remainder of the allowed keyboard 
responses, and their effects. Response '+' puts the DSP into mode 'Increase Volume' 1 1806, 
from which it returns to 'Stop?' query 1 1801; similarly, response '-' puts DSP into mode 
'Decrease Volume' 1 1809, and thence to 'Stop?' query 1 1801. Response 'u' puts DSP into 
'Uncorrected Audio Mode' 1 1807, and thence to 'Stop?' query 1 1801. Response 'i' puts DSP 
into mode 'Adjust dL/dx Correction' 1 1808, and thence to 'Stop?' query 1 1801. Response 'o' 
puts DSP into mode 'Adjust Offset' 1 1810, and thence to 'Stop?' query 1 1801. Response 'j' puts 
DSP into mode 'Adjust dL/dx Offset' 11811, and thence to 'Stop?' query 1 1801. Response 'm' 
puts DSP into mode 'Mute On' 11812, and thence to 'Stop?' query 1 1801. Response 'k' puts 
DSP into mode 'Adjust Linear Spring' 11813, and thence to 'Stop?' query 1 1801. Response T 
puts DSP into mode 'Turn IR Filter On' 1 1814, and thence to 'Stop?' query 1 1801. Response 'n' 
puts DSP into mode 'Mute Off 11815, and thence to 'Stop?' query 1 1801. Response 'v' puts the 
DSP into mode 'Adjust Nonlinear BEMF' 1 1816, and thence to 'Stop?' query 1 1801. Response 
'd' puts DSP into mode 'Turn IR Filter Off 1 181 7, and thence to 'Stop?' query 1 1801. And 
finally, a User response 's' puts the DSP into 'Stop' mode 11818, from whence the system is 
returned to 'Stop?' query 1 1801. It should be noted that all processes within the Wait Loop and 
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Command Parser, are interruptible by ISR 1 1303. The C programming language code 
implemented in the DSP is provided in file 071 1 19.txt on the compact disk s of the computer 
program listing which is a part of this application. 

DETAILED DESCRIPTION 11: METHODS AND CIRCUITS 

Page 107, line 7, please amend the following paragraph: 

In particular^ the position-indicator state variable^ /(^)? utilized in this embodiment of 
the invention is an output voltage derived from the functional dependence of the effective 
complex coil impedance Z^{(o,x) upon coil/diaphragm position x, at some fixed supersonic 

probe frequency co . The physical effects whieh -that give rise to this functional dependence, 
along with a mathematical model developed to simulate them, in accordance with the present 
invention, are described in Details 1 and 6. This embodiment is called the method.^ m-ln this 
section we elaborate on the methods and circuits used to implement the Z^ method. 

Page 108, line 4, please amend the following paragraph: 

In the context of an audio transducer, the input signal to the voice coil will include audio 
information (program material) together with the probe tone. It is therefore necessary to separate 
the probe tone and program material in frequency, so that the probe tone measurement is not 
interfered with by the audio drive signal. The Nyquist eriteri acriterion suggests that the probe 
tone 12101 should have a fi-equency of at least twice the audio fi^equency bandwidth, to avoid 
aliasing with the program material. A probe tone having a frequency of 43kHz has been found to 
be particularly desirable. However, many other frequency values could be used. 

Page 109, line 18, please amend the following paragraph: 

FIG. 55 shows a circuit diagram of the summing circuit 12202. The audio drive signal 
12201 is provided as an input to filter 12401 which isolates the probe tone 12101 from the low 
impedance of the audio amplifier output. The filter 12401 is composed of resistive, 



provided to a capacitor 124C4, which in turn is connected to the summing point 12402. 
Capacitor 124C4 decouples the audio drive signal at the summing point 12402 from the low 
impedance output of the sine wave generator 12304. The signal at the summing point 12402 is 
provided at output terminal 12403 which is connected to an input of the potential divider circuit 




^ capacitive , and inductive elements, as indicated in FIG. 55. The probe tone 12101 is 



12203. 
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Page 1 10, line 21, please amend the following paragraph: 

FIG. 58 shows the circuit of the low pass filter 12206. The first part of the low pass filter, 
incorporating the operational amplifier 1270A1, is a standard 2nd order Butterworth low pass 
filter. The second part of the filter is an inverting amplifying stage wh4eh- that includes an 
operational amplifier 1270A2 and a variable resistance 127VR1 that produces a DC offset in the 
output signal. This offset is set to reduce the DC offset in the magnitude of the probe tone that 
has been detected. The gain of the inverting amplifying stage is set to enhance the signal 
significance when it is converted to digital form. One skilled in the art will recognize that many 
different circuit arrangements could be used for the filter, gain and offset circuit, and that the 
rather straightforward circuit shown in FIG. 58 can be modified without changing the essence of 
the design. Operational amplifiers 1270A1 and 1270A2 may be National Semiconductor part 
number LM741. 

Page 111, line 18, please amend the following paragraphs: 

The filter based method used in the- Z^(x) detection circuit 12200 and shown in FIG. 54, 

is sensitive to changes in output impedance of the audio amplifier 12303. For example, with a 
low impedance load, some types of amplifiers exhibit large crossover distortion effects, which in 
effect are a change in output impedance. This change in output impedance can cause noise in the 
Z^{x) measurement. Furthermore, in transducers driven with large currents there can be 

considerable heating effects in the coil. This produces a change in the Ohmic resistance 

whte hthat is misinterpreted by the Z^{x) detection circuit 12200 as a change in position (this is 

discussed in Detail 6 above). Someone skilled in the art would recognize that a more complex 
circuit is required to separate out these two effects for the full range of transducers, but that this 
would not materially change the invention detailed here. 

It will be apparent to those skilled in the art that the particular position-indicator state 
variable / (x) described in this section and in Detail 6, which is derived from the functional 

dependence of the effective complex coil impedance Z^(co,x) upon coil/diaphragm position x 

at some fixed supersonic probe fi-equency o) , can be used within various embodiments of a 
feedback linearization control system according to the present invention, in which the positional 
information / (jc) is used in various different ways, including but not limited to one or more of 
the control laws presented in Details 2 and 10 above. 
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Page 112, line 27, please amend the following paragraphs: 

FIG. 34 shows a schematic cross section of a typical cell phone speaker or receiver 
13100; actual three-dimensional speaker geometry is a figure of revolution about the central 
horizontal axis of symmetry (not shown). Speaker 13100 consists of a transducer and integral 
acoustic venting. A voice coil 13101 is mounted on the diaphragm 13102. Coil 13101 is 
positioned in the gap between a neodymium magnet 13103 and a magnetic base plate 13104. A 
plastic surround 13105 supports the diaphragm 13102 and a face plat e faceplate 13106. The 
surround and face plate faceplate have acoustic vents 13107 which tune the frequency response of 
the speaker 13100. The depth, indicated in FIG. 34 by Dl, is typically 2 mm. The main 
difference between this type of transducer assembly and the transducer shown in FIG. 3 is the 
single surrounding support of the relatively flat diaphragm 13102. This means that the system is 
resistant to the tilt ("canting") wtrieh -that can complicate capacitance position-sensing methods 
in other transducers as described in Detail 7. 

The preferred method of detecting the variation of capacitance with coil/diaphragm axial 
position, C^^^^^.^.^ (x) , is to place the capacitance within an oscillator circuit. Changes in 

Cparasitici^) ^o changcs in coil position are then the cause of changes in the oscillator 
frequency. A frequency-to-voltage converter is then used to yield a varying signal which is a 
function of the parasitic capacitance. The varying signal can be identified with C ^^^^^.^.^(x) in 

suitable units. Thus^ as defined^ C^ro^/zc (^) identified with the position-indicator state 

variable^ /(^).? 

FIG. 60 shows a schematic of the capacitance detector and speaker arrangement, together 
with the DSP used for correction. An analog audio signal, provided over input line 13201, is 
digitized by DSP based mixed-signal controller 13202. Mixed-signal controller 13202 is 
embodied by a AD21992 chip which includes an ADC (analog voltage to digital) converter). The 
output of the DSP based controller 13202 is connected to a standard DAC (digital to analog 
voltag e ) converter) 13203. The output of the DAC 13203 is amplified by a DC-connected audio 
amplifier 13204. The output of amplifier 13204 has a drive connection 13205a to one terminal of 
the voice coil 13101 of the speaker 13100. The magnetic base plate 13104 of the speaker 13100 
has a connection 13207 to one input of an oscillator circuit 13208 (detailed in FIG. 61). Another 
input to the oscillator circuit 13208 is connected to the drive connections 13205a and 13205b of 
the coil 13101 through blocking capacitors 13209a and 13209b, respectively. The output of 
oscillator circuit 13208 is connected to a frequency to voltage converter 13210, which converts 
the variable frequency received from the oscillator circuit 13208, and also amplifies and leveN 
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shifts the varying voltage output. The output 13404 from the frequency to voltage converter 
13210, which is a measure of C^^,^„^,^ (x) (abbreviated as Cp{x) in the Figure), and hence the 

position-indicator state variable^ /(jc), is input into the mixed signal DSP controller 13202. 
Inside DSP 13202, both the analog output voltage from 13210 and the analog input audio signal 

13201 are converted into digital signals, and combined by the DSP 13202 to yield the digital 
output 1321 1 of the DSP 13202. The purpose of the DSP functionality within the controller 

13202 is to furnish the DAC 13203 with a digital signal such that the output of DAC 13203, after 
amplification by amplifier 13204, will feed the speaker-transducer voice coil with a voltage 
signal including both the audio program and a pre-distortion calculated to cancel out a significant 
portion of the nonlinearities introduced by the transducer in the course of its normal uncorrected 
operation. This e ff e ct of th e position s e nsor analog signal 13 4 01 f e d by frequ e ncy to voltag e 
conv e rt e r 13210 into the mixed signal controller 13202, is t e rm e d feedback linearization. 

Page 1 15, line 17, please amend the following paragraphs: 

FIG. 62 shows the detailed circuitry of the frequency to voltage converter 13210. 
Frequency to voltage converter 13210 consists of two parts: a frequency to pulse converter 
circuit 13401, and a low pass filter, amplifier and level shifter circuit 13402. The frequency to 
pulse converter 13401 consists of a mono-stable multi^vibrator circuit 13407 Vi4H€h -that includes 
an industry standard multi-vibrator whieh -that mav be, for example, a 74LS123 as used in this 
embodiment. The mono-stable multi-vibrator circuit 13407 takes the square wave output signal 
13304 received from the oscillator circuit 13208, which has a constant rms value, and converts 
it to a pulse train which that is provided on line 13403. The pulse train 13403 has an rms value 
varying with frequency, which is a function of the transducer coil/core capacitance C^, which in 

turn varies with coil/diaphragm position x . The low pass filter, amplifier and level shifter circuit 
13402 converts the pulse train on line 13403 to a varying analog voltage output provided on line 
13404. This varying analog voltage on line 13404 represents the varying capacitance C^(jc) . The 

low pass filter, amplifier and level shifter circuit 13402 includes an operational amplifier 
1340A1, which receives the output signal on line 13403 and, using a gain of 10 as determined by 
resistor values, low-pass-filters and offsets the signal 13403; and operational amplifier 1340A2, 
which has a gain of unity and implements a second-order Butterworth filter. These operational 
amplifiers may be embodied, for example, as National Semiconductor part number LM741, or 
equivalent. Resistor 134VR1 is adjusted such that the coil/diaphragm equilibrium position 
produces a zero output voltage. Operational amplifier 1340A2 receives, at its input terminal 
13406, the output signal provided at output terminal 13405 of operational amplifier 1340A1, and 
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then converts that signal to a voltage which is provided on line 13404 to mixed signal DSP 
13202. 

Page 116, line 30, please amend the following paragraphs: 

The particular position-indicator state variable f{x) described in Detail 7 and in this 

section, which is derived from the internal parasitic capacitance C^^^^^.,,.^ between the transducer 

voice coil and the transducer magnetic pole structure, can be used with various embodiments of a 
feedback linearization control system according to the present invention, in which the positional 
information f{x) is used in various different ways, including but not limited to one or more of 
the control laws presented in Details 2 and 10 above. 

Page 1 19, line 6, please amend the following paragraphs: 

The IR-LED 14201 and PIN diode 14202 are epoxied side-by-side onto the transducer 
frame 14203, with both diodes pointing at a reflecting region 14204 on the transducer cone 
14205. Reflecting region 14204 should subtend a sufficient angle such that, as the transducer 
cone moves, the PIN diode 14202 detector admittance cone is always pointed within the region. 
The diodes are preferably inclined towards each other and pointed towards the axis of the 
transducer at approximately a right angle to the direction of motion, or towards the curve of the 
cone. As was noted in Detail 8 above, the PIN diode output is not completely linear with cone 
position and therefore requires calibration by comparison with a metrology system. The 
position-indicator variable, / (jc) , and the degree of its non-linearity, can be varied by changing 
the positions and orientations of the two diodes relative to each other and to the transducer cone. 
Thus, there is some variation from one implementation to another and some adjustment by trial 
and error may be necessary. 

Page 120, line 3, please amend the following paragraphs: 

The following sequence of adjustments is iterated five to seven times, until the output 
waveform 14401 is about 90% of peak A/D limit: 

(a) Increase IR-LED current by adjusting variable resistor 144VR1, and thus 
output power, until the magnitude of output signal 14401 is at the limit on one excursion; 

(b) Adjust the balance by changing variable resistor 144VR2 until there is no 
output signal at terminal 14401; 
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(c) Adjust the gain of amplifier 1440P2 using variable resistor 144VR3 for 
desired peak-to-peak voltage corresponding to full motion of the transducer cone; 

(d) Turn off the coil current, readjust the balance using variable resistor 
144VR2, and zero the signal 14401 when the transducer cone is at the equilibrium point. 

DETAILED DESCRIPTION 14: IR RESULTS 

Page 120, line 24, please amend the following paragraphs: 

Audio distortions were measured, both with and without the correction, by applying an 
industry-standard two-tone SMPTE test, with audio input consisting (instead of the CD player) 
of a 60 Hz tone in conjunction with a 3 kHz tone. All four corrections described in Detail 2 
were applied by the DSP-based controller: transducer correction (spring correction S and motor 
factor correction B ), the BEMF correction, and the position dependent inductive correction. 

FIG. 65 shows a portion near 3 kHz of the FFT power spectrum distribution of the SPL 
(sound pressure level) wave-pattern picked up by a microphone in the acoustic near-field; both 
corrected spectra which is indicated by reference character 1521, and uncorrected spectra which 
is indicated by reference character 1522 are shown, and it is clearly seen that the powers in the 
60 /fe -spaced lattice of intermodulation frequency peaks, are significantly reduced when the 
correction is applied. FIG. 66 shows the low-frequency portion of the same power spectrum 
distribution, showing multiple harmonics of the 60 Hz tone; again, spectra are depicted both 
with and without correction, and again, significant reduction in the magnitude of the harmonic 
distortion peaks can be seen. 
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